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Although  the  work  on  this  program  officially  started  in  March, 
i960,  little  effort  could  be  devoted  to  it  until  about  July  1,  because  of 
staffing  difficulties.  For  this  reason,  the  first  quarterly  progress 
report  was  delayed  until  the  present  time.  This  report  follows  the 
outline  set  forth  in  the  statement  of  work  and  consists  of  three  main 
parts: 

1.  Determination  of  the  Physical  Properties  of  Hydrogen, 

2.  Cryogenic  Instrumentation,  and 

3.  Design  Principles  and  Materials  Utilization. 


1.  Physical  Properties  of  Fluid  Hydrogen 


1.  1  P-V>T,  specific  heat,  enthalpy,  entropy,  heat  of 
vaporization. 

Dr.  R.  D.  Goodwin  (Project  Leader),  D.  E.  Diller, 

H.  M.  Roder,  Dr.  L.  A.  Weber,  Dr.  B.  A.  Younglove 


CEL  started  work  on  the  tasks  in  this  section  on  March  1, 

1958.  R.  B.  Scott  and  R.  D.  Goodwin  carried  out  the  early  planning 
of  the  experimental  program.  The  setting  up  of  the  laboratory  was 
done  by  R.  D.  Goodwin  with  help  from  various  assistants,  not  all  of 
whom  are  listed  above.  The  immediate  effect  of  NASA  support  has 
been  to  increase  the  staff  from  three  to  five  scientists.  Although 
recruiting  efforts  began  early,  the  two  additional  scientists  could  not 
be  obtained  until  September,  I960.  Thus  NASA  support  of  this  project 
effectively  began  only  one  month  ago.  In  order  to  convey  the  entire 
backgroimd  of  this  project,  a  review  from  its  beginnings  will  now  be 
presented. 


First,  we  direct  attention  to  a  survey  of  the  published  data  on 
the  thermal  properties  of  fluid  hydrogen  at  low  temperatures  which  is 
of  value  as  a  guide  to  the  areas  where  more  experimental  research 
is  needed.  Because  of  its  length,  this  survey  which  has  been  extracted 
from  an  earlier  report,  has  been  placed  in  an  Appendix.  It  should  be 
noted  that  the  extensive  P-V-T  data  referred  to  are  for  normal  hydrogen. 
The  extent  of  coverage  is  indicated  in  Figure  1.  A  more  detailed 


correlation  of  the  P-V-T  data  for  normal  hydrogen  has  been  published 
by  Stewart  and  Johnson^  of  this  laboratory.  The  present  investigation 


R.  B.  Stewart  and  V.  J.  Johnson,  Advances  in  Cryogenic  Engineering, 
Vol.  5,  pp.  557-565,  1959 
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is  aimed  at  para-hydrogen,  the  form  that  is  produced  for  large- 
scale  cryogenic  uses. 

The  experimental  apparatus  at  CEL  consists  of  two  major 
parts:  (1)  the  P-V-T  apparatus  and  (2)  the  calorimeter  for  measur¬ 
ing  specific  heat  at  constant  volume  and  the  latent  heat  of  vaporization. 
Of  these,  the  former  is  completed  and  is  producing  data,  while  the 
latter  is  much  less  far  along.  The  progress  as  of  a  year  ago  and  as  of 
now  is  indicated  in  Figures  2  and  3,  respectively.  The  calorimeter 
task  has  been  inactive  in  the  last  year.  The  NASA  support  will  permit 
resumption  of  this  task  and  carrying  on  the  construction  of  apparatus 
for  it  simultaneously  with  the  accumulation  of  P-V-T  data. 

The  experimental  tasks  are  being  carried  out  with  what  are 
conceived  to  be  the  best  available  techniques  consistent  with  the  time 
schedule  that  has  been  adopted.  We  can  hardly  hope  to  convey  in 
these  reports  the  carefvtl  attention  to  detail  that  is  required  in  such 
accurate  experimentation  and  must  be  content  with  brief  summaries 
of  the  principal  features. 

The  P-V-T  apparatus  consists  of  a  small  pressure  vessel 
mounted  in  a  cryostat  and  equipped  with  precise  means  for  regulating 
and  measuring  its  temperature.  It  connects  via  a  capillary  amd  a 
pressure-balancing  diaphragm  cell  with  an  accurate  dead-weight 
rotating-piston  pressure  gage.  Temperatures  down  to  20* K  and 
pressures  up  to  5000  psi  can  be  attained  and  accurately  measured. 
Resiilts  of  one  "pseudo-isochore"  for  normal  hydrogen  are  shown  in 
Figure  4.  The  resulting  compressibility  factors  agree,  on  the  average, 
within  0.  1%  with  the  best  published  values.  Following  this  experiment, 
which  was  for  the  purpose  of  testing  the  apparatus,  sixteen  pseudo¬ 
isochores  of  para-hydrogen  have  been  obtained  at  molal  volumes  from 
28  to  40  cm^.  Further  measurements  will  be  made  at  both  higher  and 
lower  densities.  In  looking  ahead  to  the  problem  of  how  to  correlate 
a  wide  range  of  P-V-T  data,  Goodwin  has  developed  am  approximate 
empirical  equation  of  state.  ^ 


R.  D.  Goodwin,  Advances  in  Cryogenic  Engineering,  Vol.  6, 
Paper  G-4,  I960  (in  press) 
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X 
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Cryostat  electric  wiring 
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A  preliminary  design  of  the  calorimetric  apparatus  was  made 
at  an  early  stage.  This  was  desirable  because  it  was  intended  to  use 
some  of  the  same  temperature  control  and  measurement  equipment 
and  sample -handling  apparatus  for  both  the  P-V-T  and  calorimetric 
experiments.  Thus  some  degree  of  coordination  in  the  early  planning 
was  required.  Steel  hemispheres  for  the  pressure  vessel  were 
acquired,  but  construction  of  the  calorimeter  was  then  set  aside  in 
order  to  progress  with  the  P-V-T  apparatus.  In  the  past  month,  the 
calorimeter  has  again  been  taken  up  and  will  now  be  completed.  The 
pressure  vessel  is  being  assembled,  and  principles  of  thermal  shield¬ 
ing  are  being  examined.  An  early  schematic  diagram  of  the  calorimeter 
is  shown  in  Figure  5  and  will  serve  to  suggest  the  general  principles 
to  be  used. 


1.  2  Thermal  Conductivity 

W.  J.  Hall  and  R.  L.  Powell 

This  task  also  was  initiated  in  March,  1958.  A  guarded  flat- 
plate  configuration  was  adopted  for  the  experimental  cell.  As  in 
Section  1.  1,  temperatures  down  to  20“K  and  pressures  up  to  5000  psi 
can  be  attained.  Figures  6  and  7  show  the  apparatus  schematically, 
while  Figures  8,  9,  and  10  are  photographs  of  the  pressure  bomb, 
the  gas  handling  and  purification  apparatus,  and  the  electrical  instru¬ 
mentation  console,  respectively. 


The  thermal  conductivity,  K,  is  defined  by  the  Fourier  equation: 


dT 


when  W  is  the  heat  flow,  A  is  the  cross-sectional  area,  and  dT/dX  is 
the  temperature  gradient  in  the  direction  of  heat  flow.  Since  the 
measurements  will  be  made  at  equilibrium,  and  since  the  property 
being  measured  is  macroscopic  instead  of  microscopic,  the  thermal 
conductivity  is  approximated  by  the  equation: 

W  (AX) 

A  (AT) 

where  W  is  now  the  power  put  into  the  hot  plate  (See  Figure  6).  This 
is  determined  by  measuring  the  current  flowing  in  the  hot  plate  heater 
and  the  voltage  impressed  upon  it.  The  cross  sectional  area.  A,  and 
the  distance  between  the  hot  and  cold  plates,  AX,  are  determined  at 
room  temperature  and  corrected  for  temperature  and  pressure. 
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FIGURE  7 
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The  temperature  difference  between  the  hot  plate  and  the  cold  plate, 

AT,  is  measured  both  by  platinum  resistance  thermometers  and  by 
difference  thermocouples. 

Heat  flow  from  the  hot  plate  will  be  nearly  linear  into  the  cold 
plate  since  the  hot  plate  is  surrounded  by  a  guard  which  will  be  held 
at  the  same  temperature  as  the  hot  plate  to  within  0.  001*K  by  a 
thermopile  composed  of  eight  difference  thermocouples.  The  auxiliary 
guard  is  needed  to  prevent  excessive  cooling  of  the  guard  by  cold  gas 
flowing  into  the  pressure  cell  from  above. 

The  hot  plate  is  supported  in  three  places  by  quartz  spacers 
forced  into  it  by  stainless  steel  screws.  After  the  hot  plate  and  guard 
were  fastened  together  they  were  machined  smooth  on  the  bottom. 
Fused  quartz  spacers  placed  under  the  guard  support  both  the  guard 
and  the  hot  plate  above  the  cold  plate.  Several  values  of  AX,  e.  g. 

0.  25,  0.  50,  1.00  mm.  are  available. 

The  temperature  of  the  hot  plate  and  cold  plate  are  measured 
by  platinum  thermometers  imbedded  in  them  and  open  to  the  pressure 
in  the  cell.  Since  the  difference  in  temperature  is  more  easily 
stabilized  from  thermocouple  output,  thermocouples  have  been 
installed  in  parallel  with  the  platinum  thermometers.  The  absolute 
temperature  is  measured  by  a  sealed  platinum  thermometer  imbedded 
in  the  wall  of  the  pressure  cell. 

It  is  possible,  in  fact,  necessary,  to  determine  the  pressure 
coefficients  of  resistance  for  the  unsealed  platinum  thermometers  and 
the  thermocouples  used.  This  can  be  accomplished  with  this  arrange¬ 
ment  since  the  cold  plate  platinum  thermometer  and  the  calibrated, 
sealed  platinum  thermometer  in  the  cell  wall  are  well-connected 
thermally. 

To  keep  the  measuring  components  as  isothermal  as  possible, 
they  have  been  fabricated  from  electrolytic  tough  pitch  (99.  9  per  cent 
pure)  copper.  It  was  necessary  to  make  the  pressure  cell  of  beryllium 
copper  to  contain  the  pressure  even  though  its  thermal  conductivity  is 
undesirably  low.  To  improve  the  heat  flow  pattern,  a  cooling  block  is 
attached  to  the  bottom  of  the  pressure  cell  (See  Figure  6).  Heaters 
and  difference  thermocouples  were  put  in  various  places  on  the 
pressure  cell  to  make  it  as  isothermal  as  possible. 
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The  measuring  components  and  associated  pressure  container 
are  supported  from  above  by  a  stainless  steel  tube  through  which  the 
electrical  wiring  enters  the  measuring  cell.  This  tube  thermally 
isolates  the  pressure  cell  from  the  boiling  cryogen.  By  controlling 
the  amount  of  boil-off  gas  passed  through  the  cooling  coils  and  the 
power  input  to  the  external  heaters,  the  measuring  cell  can  be  raised 
to  desired  temperatures  above  the  boiling  point  of  the  cryogen.  The 
large  heat  capacity  of  the  cell  facilitates  temperature  stabilization. 

The  foregoing  apparatus  was  nearly  all  completed  prior  to 
initiation  of  NASA  support.  During  the  period  covered  by  this  report, 
calibration  of  the  platinum  resistance  thermometers  within  the  gas 
space  of  the  cryostat  has  proceeded  smoothly  and  regularly.  Measure¬ 
ments  have  been  made  at  five  well -separated  temperatures  and  at 
pressures  of  0.  1,  1,  10,  50,  100  and  250  atmospheres.  As  expected, 
the  resistances  of  the  thermometers  vary  linearly  with  temperature, 
so  that  calibration  and  interpolation  calctilations  are  both  simple  and 
accurate.  The  pressure  changes  during  an  experimental  run  are 
vindetectable  for  the  lower  pressures  and  only  about  one  part  in  a 
thousand  for  the  highest  pressures.  The  temperature  drifts  of  the  ^ 
cold  and  hot  plates  within  the  cryostat  are  held  to  about  1  to  3  x  lO"  *K 
for  each  run.  A  refinement  in  measuring  techniques  has  pushed  down 
the  uncertainties  in  measurement  of  temperature  differences  to  about 
the  same  level.  A  servo-control  for  the  power  supply  to  the  main 
block  heater  has  brought  about  both  a  more  rapid  approach  to  steady 
state  conditions  and  also  smaller  cyclic  variations  in  the  block 
temperature. 

The  calibration  runs  have  had  the  special  additional  benefit  of 
revealing  fine  details  of  the  temperature  distributions  within  the  system. 
It  was  found  recently  that  there  were  gradients  of  up  to  0.  01*K  between 
the  cold  plate  and  the  block  and  that  these  gradients  were  pressure 
dependent.  The  expected  limit  of  accuracy  for  the  system  as  a  whole 
is  also  about  0.  01*K.  However,  the  gradients  mentioned  would  have 
no  effect,  in  the  first  approximation,  on  actual  thermal  conductivity 
readings.  Nevertheless,  the  system  is  being  changed  slightly  in  order 
to  reduce  minor  fluctuations  in  the  later  measurements. 

The  gradients  were  caused  by  the  cooling  of  the  outer  guard  and 
the  subsequent  cooling  of  the  inner  guard,  hot  plate,  and  cold  plate. 

The  cooling  of  the  outer  guard  was  caused  by  conduction  up  the  copper 
wires  and  by  gas  conduction  to  the  upper  plate  of  the  large  pressure 
block.  The  gradients  observed  were  pressure  dependent  because  the 
gas  conduction  from  the  outer  guard  to  the  inner  parts  was  itself 
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pressure  dependent.  The  new  additions  to  the  cryostat  are: 

(1)  a  heater  on  the  outer  guard  to  balance  the  heat  conduction  up  the 
wires;  and  (2)  a  heater  on  the  top  plate  of  the  block  to  bring  it  into 
thermal  equilibrium  with  the  main  part  of  the  block.  The  necessary 
auxiliary  power  supplies  and  controllers  are  being  placed  on  the 
instrument  panels. 

Following  completion  of  the  thermometer  calibrations,  a  brief 
series  of  measurements  will  be  made  on  air  at  low  pressures  in 
connection  with  another  project.  The  apparatus  will  then  be  devoted 
to  making  the  more  extensive  high-pressure  measurements  on  para- 
hydrogen. 


1.  3  Dielectric  Constant 


Experimental  measurements  of  dielectric  constant  are  being 
deferred  for  the  present  for  the  following  reasons: 

a.  Examination  of  the  published  data  for  normal  hydrogen 
shows  that  the  Clausius -Mosotti  equation  fits  the  best  data  within 
about  1%  in  (e-  1)  over  a  very  great  range  of  conditions  from  low 
pressure  gas  at  ambient  temperature  to  saturated  liquid.  This  indicates 
that  only  a  few  experimental  measurements  will  need  to  be  made  in 
order  to  establish  the  pattern  of  the  deviation  (if  any)  from  the 
equation.  Interpolation  throughout  the  P-T-p  region  of  interest  can 
then  be  made  using  the  equation  rather  than  by  direct  measurement. 
Since  the  experiment  will  be  a  relatively  simple  one,  it  does  not  have 
to  be  started  at  once.  Also,  the  correlation  of  the  results  would  have 
to  wait  anyway  on  the  availability  of  accurate  data  from  the  P-V-T 
project. 


b.  At  a  later  date,  it  may  be  possible  to  use  a  pressure  cell, 
cryostat,  and  gas  handling  equipment  developed  for  one  of  the  other 
tasks,  such  as  the  one  on  thermal  conductivity,  and  thus  avoid  duplica¬ 
ting  this  expensive  equipment. 

1.4  Sonic  Velocity 

We  have  attempted  vigorously  to  recruit  a  man  with  experience 
in  ultrasonics  for  this  task,  but  have  not  been  successful  as  yet.  Some 
good  prospects  were  lost  because  the  recent  federal  pay  raise  and 
liberalization  of  Civil  Service  hiring  procedures  for  top-quality  Ph.  D.  's 
came  a  little  too  late  relative  to  the  period,  February-May,  when  most 
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recruiting  is  accomplished.  Recruiting  efforts  will  continue,  of 
course,  with  the  aim  of  providing  a  top-quality  project  leader. 
Consultative  services  in  ultrasonics  are  already  available  within 


NBS. 


2.  Cryogenic  Instrumentation 

J.  Macinko,  P.  Smelser,  R.  C.  Muhlenhaupt, 

C.  E.  Miller,  R.  B.  Jacobs 

The  first  phase  of  the  Cryogenic  Instrumentation  project  is 
the  determination  of  the  status  of  instrumentation  in  the  areas  of 
pressure,  temperature,  liquid  level  and  flow  rate  measurement.  A 
survey  of  technical  publications  revealed  that  very  little  had  been 
published  with  regard  to  response  and  performance  characteristics; 
the  btilk  of  the  literature  was  descriptive  in  nature  but  in  many  cases 
the  theory  of  operation  was  also  covered. 

Letters  of  inquiry  were  sent  to  industries  and  laboratories 
known  to  be  working  in  cryogenics  requesting  information  on  the  use 
of  instruments  in  their  facilities.  The  answers  received  were  helpful 
in  providing  direction  for  future  test  programs.  However,  many 
inquiries  remain  unanswered  and  many  answers  were  uninformative; 
the  overall  result  was  disappointing . 

Several  reports  of  the  results  of  test  and  evaluation  programs 
were  also  studied.  Some  of  these  reports  are  excellent  in  content 
and  organization.  This  study  indicates  that  broader  coverage  of  the 
field  is  definitely  required.  In  addition,  testing  in  the  cryogenic  region 
has  been  very  limited.  Manufacturer's  literature  and  performance 
specification  sheets  provide  additional  information,  but  in  most  cases 
the  manufacturers  have  not  provided  performance  characteristics  in 
temperature  and  pressure  regions  of  importance  to  the  cryogenic 
industry. 

Due  to  the  above  factors,  the  first  phase  of  the  instrumentation 
survey  contains  very  little  evaluation  based  on  controlled  tests  or  user 
application.  However,  a  compilation  of  the  various  methods  of  measure¬ 
ment  of  the  four  parameters  under  study,  the  physical  principles 
involved,  and  a  partial  list  of  manufacturers  of  the  instrtiments  under 
investigation  may  be  of  immediate  value  to  engineers,  scientists  and 
technicians  engaged  in  cryogenics.  This  compilation  is  nearing 
completion  and  has  been  set  up  using  an  orderly  index  which  permits 
additions  under  all  classifications  and  subheadings.  Instrvunents  have 
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been  evaluated  wherever  experience  or  data  permitted.  Additions 
and  revisions  will  be  made  as  test  results  and  manufacturers' 
literature  are  received. 

A  test  program  is  planned  for  all  phases  of  the  study  and  will 
be  entered  into  progressively.  An  apparatus  has  been  designed  and 
built  for  testing  pressure  transducers  (Figure  11)  and  will  soon  be 
ready  for  its  initial  cooldown  and  performance  evaluation.  An  idea 
for  a  densitometer  is  being  studied  and  a  model  has  been  built  and 
will  be  tested  within  the  next  month  (Figure  12).  Several  ideas  have 
been  investigated  for  determining  response  time  of  temperature 
sensors  and  final  plans  of  a  test  apparatus  will  be  drawn  up  shortly. 

No  additional  test  plans  have  been  formulated  to  date,  but  additional 
testing  will  be  planned  as  soon  as  our  studies  indicate  which  parame¬ 
ters  should  be  measured,  the  conditions  under  which  measurements 
shoiild  be  made,  and  the  degree  of  accuracy  required  in  testing. 

2.  1  Forced  Vibration  Densitometer  Study 

An  experimental  program  is  being  conducted  to  determine  the 
feasibility  of  using  a  vibrational  system  as  a  means  of  determining 
flvdd  densities.  The  underlying  idea  is  that  the  mass  of  a  vibrating 
system  is  a  primary  factor  in  determining  the  dynaucnic  characteristics 
of  the  system.  If  certain  characteristics  of  the  system  are  appreciably 
affected  by  the  density  of  the  fluid  flowing  through  the  system,  a 
properly  instrumented  system  will  provide  a  means  for  measuring 
density. 

The  choice  of  the  system  has  depended  on  many  things,  among 
which,  and  perhaps  most  important,  is  the  reliability  and  accuracy  of 
associated  mechanical  and  electronic  instruments. 

Figure  12  illustrates  the  system  that  is  presently  being  investi¬ 
gated.  The  flow  passage  is  supported  by  bellows  with  an  effective 
transverse  spring  constant  kg.  The  passage  is  being  driven  with  a 
scotch  yoke  mechanism  which  has  been  designed  to  transmit  a  sinusoi¬ 
dal  motion  with  an  amplitude  sxifficiently  large  to  influence  the 
fluid.  The  motion  of  the  drive  mechanism  is  transmitted  through  a 
force  gage  to  the  passage.  A  force  transducer  with  an  extremely  high 
spring  constant  has  been  selected.  By  mcd^ing  the  spring  constant 
large  enough,  we  can  make  the  assumption  that 


^(t)  ■  ^(f) 
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where  X.  .  is  the  transverse  motion  of  the  passage  and  X.,.  the  forcing 
function.  ^Writing  the  equation  of  motion  of  the  system  we  nave 

MX  +  CX  +  K  X  =  F 
B 

where  M  is  the  total  mass  of  the  system  (M  =  m+pV),  and  m  is  the 
effective  mass  of  all  the  vibrating  mass  except  the  fluid,  V  is  the 
volume  of  the  passage,  and  p  the  density  of  the  fluid.  Some  damping 
is  inherently  involved  in  this  type  of  system  and  is  shown  as  CX. 

It  can  be  shown  from  the  equation  of  motion  that  by  measuring 
the  force  being  transmitted  to  the  passage  at  the  point  of  maximum 
acceleration  that  the  density  parameter  can  be  attained  and  is 
related  to  the  force  by  the  following  equation 

p  =  A  +  BF 

where  A  and  B  are  constants  that  can  be  determined  experimentally. 

3,  Cryogenic  Design  Principles  and  Materials  Utilization 
D.  B.  Chelton  and  B,  W.  Birmingham 

The  statement  of  work  calls  for  a  survey  to  be  initiated  by 
NBS  in  the  broad  area  of  cryogenic  design  principles  and  materials 
utilization  to  provide  support  for  the  Centaur  and  Rover  program.  It 
has  been  difficult  to  know  exactly  how  such  a  task  should  be  conducted. 
To  get  things  started,  contact  has  been  made  with  Convair  Astronautics 
regarding  Project  Centaur  and  the  Los  Alamos  Scientific  Laboratory 
regarding  Project  Rover.  Two  visits  were  made  to  Convair  Astro¬ 
nautics  and  several  telephone  discussions  have  taken  place.  In 
addition,  personnel  from  Convair  Astronautics  have  visited  CEL  on  at 
least  two  separate  occasions.  Telephone  conversations  have  also 
taken  place  with  LASL  personnel,  and  members  of  LASL  divisions  N, 

J,  and  CMF  visited  CEL  to  discuss  cryogenic  engineering  problems 
on  Rover  and  determine  what  research  was  currently  underway  at 
CEL.  Some  progress  has  been  made  in  this  manner  both  in  learning 
about  problems  and  becoming  acquainted  with  project  personnel. 

A  detailed  description  of  some  of  the  technical  problems  that 
have  been  considered  on  Project  Centaur  follows: 

3.  1  Propellant  Tank  Insulation 

The  insulation  of  the  Centaur  propellant  tank  has  been  studied 
and  several  problems  encountered. 
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Improper  application  of  the  insulation  on  the  forward  bulkhead 
of  the  fuel  tank  can  condense  a  sizable  quantity  of  air.  This  can 
result  in  several  undesirable  conditions:  (a)  The  solid  and/or  liquid 
air  can  add  considerable  dead  weight  to  the  vehicle,  (b)  The  solid 
air  can  freeze  the  side  insulation  panels  in  place  so  they  cannot  be 
jettisoned  when  required,  (c)  The  condensation  process  can  add 
considerably  to  the  anticipated  fuel  evaporation  rate  by  increasing 
the  ''apparent"  thermal  conductivity  of  the  insxilation. 

The  method  of  application  of  pre-formed  insulation  does  not 
completely  cover  the  tank  to  assure  svifficient  vapor  barrier  to  pre¬ 
vent  condensation  of  air.  Also  certain  areas,  dictated  by  geometry, 
may  prevent  application  of  pre-formed  insulation.  A  suitable  bonding 
agent  may  also  represent  a  formidable  problem  area. 

The  spray  foam  method  of  applying  insulation  has  been  devel¬ 
oped  by  Convair  Astronautics  and  other  investigators  to  a  point  where 
serious  consideration  should  be  given  to  such  application  on  the 
forward  bulkhead.  This  technique  will  form  a  much  better  vapor 
barrier  especially  over  the  parts  having  diffictilt  geometrical  configu¬ 
rations.  In  the  past,  great  difficulty  has  been  experienced  with  foams 
breaking  away  from  the  surface  and  cracking  when  thermal  cycled. 
However,  it  appears,  from  the  limited  experimental  data  available 
that  Convair  Astronautics  may  have  solved  the  adhesion  and  cracking 
problems  with  their  present  spraying  technique.  These  data  are 
available  on  small  tanks  only,  and  initial  tests  showed  good  adhesion 
for  a  single  cool-down.  However,  thermal  cycling  experiments,  done 
later  at  the  suggestion  of  CEL>,  induced  failures.  Proper  surface 
cleaning  techniques  recommended  by  CEL  resulted  in  greater  success. 
Thus,  once  again  it  is  shown  that  fabrication  control  techniques  are 
essential  when  applying  these  foams. 

CEL  has  available  a  company  to  spray  foam  samples  (in 
Boulder,  Colorado)  if  additional  data  are  required.  A  limited  number 
of  tests  could  be  conducted  if  additional  information  is  needed  to  back 
up  the  work  being  done  at  Convair.  It  is  unfortunate  that  Convair  did 
not  know  of  research  that  has  been  carried  out  by  NBS,  Beechcraft 
Research  and  Development,  and  Armour  Research  Foundation.  The 
latter  organization  has  been  working  for  a  number  of  years  to  perfect 
such  insulation  methods. 
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Although  the  insulation  program  at  Convair  Astronautics  is 
meeting  with  a  good  degree  of  success,  a  more  fundamental  approach 
to  the  problem  would  seem  advisable.  The  accelerated  time  schedule 
present  in  their  program  is  well  recognized.  However,  a  parallel 
program  could  be  established  to  more  fully  understand  the  pertinent 
considerations  of  the  application.  For  instance,  it  has  been  found  by 
previous  investigators  that  slight  variations  in  foam  composition, 
surface  finish  and  preparation,  method  of  application,  spray  timing 
sequences,  etc.  ,  materially  influence  the  results.  It  is  therefore 
conceivable  that  a  success  on  one  tank  would  not  necessarily  insure 
success  on  all  tanks  unless  a  high  degree  control  is  exercised  in  both 
application  and  foam  composition.  Since  the  important  factors  of  foam 
composition  are  not  fully  understood,  the  aforementioned  quality 
control  would  seem  difficult. 

3.  2  Properties  of  Cryogenic  Materials 

General  assistance  has  been  given  to  various  groups  within  the 
Centaur  Project  on  properties  of  cryogenic  materials  --  fluids, 
structural  materials  and  insulations.  Thus  far,  the  assistance  has, 
in  general,  been  acquainting  the  personnel  with  properties  that 
already  exist  in  the  literature.  The  Cryogenic  Data  Book,  WADC  59-8 
prepared  by  the  NBS  is  in  general  use. 

a.  Fluid  Properties 

In  addition  to  acquainting  personnel  with  existing  data  from 
literature,  several  problems  were  studied  in  connection  with  the 
properties  of  hydrogen  in  its  various  ortho-para  modifications.  An 
example  is  the  study  of  several  existing  temperature -entropy  diagrams 
for  hydrogen  to  determine  which  is  the  most  appropriate.  By  basic 
thermodynamic  reasoning,  it  was  possible  to  disprove  the  validity  of 
data  in  certain  areas  of  an  often  used  diagram.  Information  was  also 
supplied  on  the  gas-liquid  phase  equilibrium  for  the  system  helium  - 
hydrogen.  (See  Section  3.  6,  Thermodynamic  Charts  for  Hydrogen. ) 

Much  information  is  lacking  in  the  low  temperature  compressed 
fluid  region  for  hydrogen.  These  properties  are  presently  being 
investigated  at  CFL  and  will  be  made  available  as  soon  as  possible. 

b.  Material  Properties 

It  is  conceivable  that  certain  mechanical,  transport  and 
expansivity  properties  of  structural  materials  applicable  to  the  Centaur 
Project  have  not  yet  been  determined.  The  test  facilities  at  CEL>  may 
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be  useful  in  assisting  the  Convair  group  in  such  determinations. 
However,  the  present  heavy  commitments  of  the  CEL  mechanical 
properties  project  will  undoubtedly  limit  extensive  testing  in  this 
area.  The  assistance  of  CEL  staff  has  already  been  used  by  Convair 
in  the  establishment  of  such  a  testing  facility.  It  is  our  understanding 
that  this  facility  is  quite  active  and  producing  useful  results. 

c.  Bearings 

A  proposed  centrifugal  zero-gravity  vent  device  requires 
bearings  which  operate  surrounded  by  low-temperature  gas  phase  fuel 
since  liquid  cooling  cannot  be  assured.  Conventional  lubrication  is, 
of  course,  impossible.  A  high  degree  of  reliability  is  essential. 
Several  types  of  bearings  may  be  practical  for  such  application. 

Three  such  bearing  types  are  being  explored  on  other  programs  at 
CEL.  Gas  lubricated  and  magnetic  bearings  are  presently  in  the 
early  stage  of  development.  Their  development  cannot  be  assured  in 
time  to  be  of  use  to  the  present  Centaur  Project,  but  may  be  appli¬ 
cable  in  later  phases. 

A  test  program  to  evaluate  ball  bearing  materials  at  cryogenic 
temperatures  has  been  underway  for  several  years  at  CEL,  However, 
the  program  thus  far  has  mainly  investigated  bearings  submerged  in 
cryogenic  liquids.  To  determine  the  feasibility  of  using  badl  bearings 
operating  in  cryogenic  gases  (particularly  hydrogen)  for  the  centrifugal 
vent  devices,  preliminary  experiments  were  performed  in  a  CEL 
bearing  tester.  As  a  first  attempt,  the  most  promising  bearing 
(determined  from  experiments  with  the  bearings  submerged  in  liquid), 
was  operated  in  a  cold  nitrogen  (100“K)  gas  environment.  These 
bearings  were  constructed  of  440  C  stainless  steel  balls  and  races  with 
Fulon  A  (reinforced  polytetrafluoroethylene  resin  which  contains  a 
highly  inert,  inorganic  silicate  base  filler)  separators.  Although  the 
tests  were  not  conclusive  since  only  two  such  bearings  were  operated, 
10  hours  of  operation  were  obtained  at  100®K  under  an  essentially  no- 
load  condition.  At  the  end  of  this  period  there  was  no  noticeable 
increase  in  operating  torque.  On  close  inspection,  the  bearings  were 
seen  to  be  slightly  roughened.  However,  it  was  concluded  that  the 
operation  of  ball  bearings  in  a  cold  gas  environment  may  be  possible. 

If  it  is  considered  advisable,  future  tests  can  be  made.  In 
such  event,  several  preliminary  ’’run-in’’  phases  are  planned  that 
should  greatly  increase  the  bearings'  useful  life.  Convair  has  been 
made  aware  of  CEL  willingness  to  perform  these  future  tests.  Since 
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the  bearing  apparatus  is  part  of  another  project,  the  tests  would  have 
to  be  coordinated  with  the  existing  program,  however,  no  particular 
delays  are  foreseen. 

d.  Adhesives 

The  use  of  epoxy  resins  for  cryogenic  applications  is  di8> 
cussed  elsewhere  in  the  report.  Their  application  at  CEL,  however, 
has  provided  sufficient  knowledge  to  include  their  use  in  many 
instances.  An  example  is  the  observation  windows  for  the  zero- 
gravity  test  container. 

3.  3  Zero-Gravity  Test  Container 

Several  discussions  were  held  both  at  Convair  Astronautics 
and  in  Boulder  concerning  the  test  container  for  the  study  of  fluid 
behavior  in  or  near  zero-gravity.  Other  than  general  cryogenic 
design  principles,  assistance  was  given  in  two  particular  areas  -- 
optics  and  seals. 

Design  criteria  for  the  successful  design  of  observation  win¬ 
dows  was  established  by  the  Bureau.  The  success  of  a  window  was 
determined  by  good  structural  integrity  and  by  the  ability  to  maintain 
glass  temperatures  above  the  dew  point.  The  basic  difference  between 
the  present  design  and  other  observation  windows  used  elsewhere  is 
reflected  in  the  need  for  having  a  support  tube  between  the  two  windows. 
Heretofore,  it  was  general  practice  at  Convair  to  have  a  helium  purge 
on  the  warmer  window  to  prevent  accumulation  of  frost.  For  the 
present  application  it  was  considered  impractical  to  have  such  a 
purge.  Actual  calculations  necessary  to  insure  a  proper  design 
necessitated  several  assumptions.  To  substantiate  the  design  and 
the  assumptions  made,  laboratory  tests  were  conducted  at  CEL. 

Three  successful  methods  of  obtaining  an  observation  window  were 
developed.  The  arran^,ement  most  applicable  for  the  subject  test 
container  consisted  of  tapered  pyrex  glass  windows  sealed  to  the 
metal  support  tube  with  epoxy  resin.  A  thin  application  of  Armstrong 
A-4  epoxy  with  Activator  E  was  found  to  work  best.  The  configuration 
used  resulted  in  a  glass  temperature  only  a  few  degrees  below  the 
environment. 

Cryogenic  seals  developed  at  this  laboratory  were  directly 
applicable  to  the  test  container.  However,  because  of  the  advanced 
stage  of  the  design,  it  seemed  impractical  to  make  use  of  these 
developments.  The  seals  used  are  not  suitable  for  an  absolute  seal. 
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but  should  be  suitable  to  prevent  major  leakage  of  hydrogen  gas. 
Although  it  should  be  noted  that  one  hydrogen  fire  developed  on  a 
Point  Loma  test  stand  because  of  a  faulty  seal. 

The  low  refractive  index  of  liquid  hydrogen  calls  for  rather 
specialized  illumination  and  photographic  systems.  For  saturated 
liquid  hydrogen  at  one  atmosphere,  the  refractive  index  is  1.  1097  for 
a  wavelength  of  0.  5790  p.  This  gives  a  large  forward  scattering  of 
light.  The  difficulties  of  right  angle  illumination  (with  respect  to  the 
camera)  may  be  illustrated  from  the  fact  that  the  scattering  of  light 
from  a  hydrogen  bubble  is  10,000  times  greater  at  2*  from  the 
direction  of  illumination  than  at  90*. 

Direct  or  right  angle  illumination  of  certain  test  containers, 
however,  often  results  in  satisfactory  photographs.  It  is  felt  that 
this  is  due  primarily  to  multiple  reflections  from  the  highly  reflective 
tank  walls.  Under  certain  conditions  (dull  tank  walls)  the  events 
within  a  vessel  may  be  obscured  by  using  such  illumination. 

3.4  Zero-Gravity  Vent  Devices 

Vent  devices  for  the  fuel  tank  under  conditions  of  zero-gravity 
have  been  considered  in  some  detail.  A  great  number  of  venting 
methods  exist,  although  these  may  basically  be  classified  as  centri¬ 
fugal  and  thermodynamic.  Convair  is  designing  and  fabricating  a 
centrifugal  device  and  is  also  considering  thermodynamic  devices. 

In  addition  to  this  approach,  an  outside  contractor  will  be  (or  has  been) 
selected  to  construct  another  centrifugal  device  of  different  fundamen¬ 
tal  design. 

The  NBS  has  inspected  the  Convair  centrifugal  device  from  a 
cryogenic  aspect.  The  foreseeable  problem  areas  have  been  mention¬ 
ed.  These  problems  are  primarily  shaft  seals  and  reliable  bearings. 
Brief  laboratory  experiments  (described  elsewhere  in  this  report) 
have  been  performed  to  study  the  feasibility  of  ball  bearings  in  such 
an  application.  Although  shaft  seals  by  others  have  performed  with  a 
good  degree  of  success,  there  is  reason  to  doubt  high  reliability.  The 
''conventional "  shaft  seal  used  for  cryogenic  applications  (booster 
pumps,  etc. )  are  usually  composed  of  graphite  riding  on  a  chrome 
plated  stainless  steel  disc.  The  wear  and  deterioration  of  graphite 
over  relatively  long  periods  of  time,  when  not  lubricated,  has  been 
observed.  A  more  extensive  investigation  to  substantiate  reliability 
may  be  beneficial. 
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3,  5  Safety 

Because  of  the  need  to  meet  definite  schedules  on  a  program, 
proper  attention  is  often  not  given  to  safety.  The  safety  procedures 
at  Convair  Astronautics  -  Centaur  Project,  and  LASL  -  Rover  Project, 
have  not  been  reviewed  in  detail  by  the  CEL.  staff.  In  a  cursory 
inspection  of  the  Convair  facilities,  no  unsafe  practices  were  observed. 
It  would  seem  advisable  however,  to  review  the  test  facilities  and 
procedures  in  greater  detail,  observe  a  tanking  operation  at  Point 
Loma  and  eventually  review  procedures  and  operations  at  other  test 
facilities  planned  by  Convair.  LASL>  personnel  have  requested  an 
independent  inspection  of  their  Rover  facilities,  to  help  guarantee 
that  all  conceivable  safety  precautions  have  been  taken. 

The  NBS  staff  has  been  called  upon  by  many  agencies  of 
Government  and  Industry  to  give  advice  on  both  new  and  existing 
hydrogen  facilities.  Our  staff  is  well  acquainted  with  safety  studies 
made  by  others  and  safety  procedures  taken  by  many  laboratories  in 
addition  to  our  own  intensive  program  to  assure  safe  operations  at 
CELi  facilities. 


3.  6  Charts  of  Thermodynamic  Properties  of  Hydrogen 
Edmund  H.  Brown 

Although  the  number  of  charts  available  on  the  thermodynamic 
properties  of  hydrogen  is  not  large,  the  design  engineer  can  obtain 
several  of  them  which  are  at  least  partially  contradictory  and  for 
which  he  usually  has  neither  the  time  nor  additional  information 
required  for  a  critical  comparison.  The  purpose  of  this  section  is  to 
make  some  brief  comments  on  both  the  published  and  unpublished 
charts  presently  being  used  or  circulated  in  the  cryogenic  industry 
with  a  view  to  assisting  the  engineer  who  is  faced  with  such  a  choice. 

Published  charts  for  normal  hydrogen  include  the  three 
temperature-entropy  diagrams  in: 

1.  Keesom  and  Houthoff.  Leiden  Comm.  Suppl.  65d  (1928) 

2.  Woolley.  Scotty  Brickwedde,  NBS  J.  Research  (1948) 

3.  W.  Koeppe.  Kaltetechnik  ^  ( 1956) 

In  addition,  the  following  unpublished  diagrams  have  had  some 
circulation  in  industry: 
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4.  Temperature-Entropy  diagram  entitled  "Hydro gen -Normal” 
NACA  File  No.  4842  (1958) 

5.  Enthalpy-Pressure  diagram  "for  parahydrogen",  no 
further  identification,  prepared  by  Dr.  H.  L.  Johnston 
(on  3  sheets  covering  the  range  0  to  200  atm.) 

A  recent  Russian  temperature-entropy  diagram  has  been  reported, 
but  whether  it  is  a  partial  revision  of  earlier  diagrams  or  represents 
new  data  is  unknown  since  copies  have  been  unobtainable. 

Most  of  the  following  comments  are  restricted  to  several  key 
points  which,  although  they  cannot  be  used  to  demonstrate  the  reliabil¬ 
ity  or  accuracy  of  a  given  chart,  can  be  used  for  the  converse  --  that 
is,  to  conclude  that  a  particular  chart  contains  sufficient  inaccuracies 
that  its  use  for  design  should  not  be  seriously  considered.  These 
points  include  the  following: 

a.  no  isochore  can  intersect  the  (Joule-Thomson)  inversion 
curve  more  than  once; 

b.  the  pressure  component  of  the  curvature  vector  of  the 
inversion  curve  in  a  In  p  -  In  T  diagram  must  be  every¬ 
where  negative; 

c.  the  low-temperature  terminal  point  of  the  inversion 
curve  must  lie  in  the  saturated  liquid  between  the  triple 
point  and  the  critical  point; 

d.  the  slope  of  the  isenthalps  in  a  T-s  diagram  must  be  every¬ 
where  negative  in  the  two-phase  (liquid-vapor  or  solid- 
vapor)  region; 

e.  both  the  temperature  and  entropy  components  of  the  curva¬ 
ture  vector  of  the  isenthalps  in  a  T-s  diagram  must  be 
everywhere  positive  in  the  two-phase  region; 

f.  the  enthalpy  of  the  saturated  vapor  must  pass  through  a 
maximum  at  some  pressure  below  the  critical  pressure. 

Some  of  the  above  conditions  are  discussed  in 

6.  E.  H.  Brown,  On  the  thermodynamic  properties  of  fluids, 
Proc.  I960  Session  Comm.  I  of  the  Int.  Inst,  of  Refr. 

(to  appear); 

the  remainder  will  be  discussed  briefly  here. 

The  early  Leiden  diagram  (1)  of  Keesom  and  Houthoff  appears 
to  satisfy  all  the  above  conditions.  Actually  there  are  two  low- 
temperature  diagrams  in  this  article;  in  the  second  (or  "corrected" 
diagram)  the  isenthalps  in  the  compressed-liquid  region  are  slightly 
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flatter  than  the  original.  In  both,  the  enthalpy  of  the  saturated  vapor 
passes  through  a  maximum  at  about  3.  S  atm.  The  temperature  at 
which  the  inversion  curve  terminates  in  the  saturated  liquid  is  approx¬ 
imately  27* K  in  the  original  diagram,  and  slightly  lower  (perhaps  as 
low  as  26* K)  in  the  "corrected"  diagram. 

The  diagram  (2)  of  Woolley,  Scott,  Brickwedde  is  generally 
similar  to  the  Leiden  diagram  (1)  and  also  appears  to  satisfy  all  the 
above  conditions.  Since  this  later  diagram  contains  considerable 
newer  data  and  is  the  result  of  a  wider  and  more  careful  correlation 
of  data,  some  preference  should  be  given  for  its  use  rather  than  (1); 
also,  it  is  more  readily  available. 

The  recent  T-s  chart  due  to  Koeppe  (3)  has  been  discussed  in 
some  detail  in  (6).  The  inversion  curve  passes  through  the  point 
55* K,  95  atm.  with  nearly  the  same  slope  as  the  500  Amagat  isochore 
which  lies  about  2*K  above  this  point;  thus,  only  a  slight  extrapolation 
to  higher  density  is  needed  to  show  that  in  this  region  conditions  (a) 
will  be  violated.  Condition  (b)  is  violated  by  all  the  points  representing 
Koeppe' a  own  data;  the  only  points  not  showing  a  positive  curvature 
are  those  from  the  (also  probably  incorrect)  correlation  by  Baehr  of 
early  J-T  measurements  of  Johnston,  Bezman,  and  Hood.  The  inver¬ 
sion  curve  appears  to  wander  off  in  the  general  direction  of  the  fluid- 
solid  phase  transformation  curve;  thus,  condition  (c)  is  probably 
violated.  Although  conditions  (d)  and  (e)  are  satisfied,  condition  (f) 
is  not  --  there  is  no  sign  of  a  maximum  in  the  enthalpy  of  the  saturated 
vapor,  at  least  to  pressures  below  one  atm.  Because  of  these  many 
inaccuracies  and  irregularities,  use  of  Koeppe's  diagram  (3)  cannot 
be  recommended. 

The  ubiquitous  NACA  temperature-entropy  chart  appears  to 
have  been  widely  circiilated  in  the  missile  industry  despite  the  fact  that 
it  is  clearly  marked  "data  not  checked"  and  contains  many  obvious 
Inaccuracies.  In  addition  to  the  210  BTU/lb.  isenthalp  being  carelessly 
drawn  and  inconsistent  with  neighboring  isenthalps,  the  following 
features  can  be  noted:  the  smallest  isenthalp  included  shows  a  clear 
maximum  at  42.  2*R  or  23.  4*K  and  the  extrapolated  terminal  point  of 
the  inversion  curve  is  at  about  35* R  or  19.  5*K,  a  temperature  far 
below  the  approximate  27 *K  indicated  by  other  reliable  charts  as  well 
as  all  the  recent  data;  both  conditions  (d)  and  (e)  are  violated  in  parts 
of  the  two-phase  region.  Because  of  the  importance  of  these  latter 
conditions  to  the  analysis  of  this  chart,  we  shall  give  a  sketch  of  their 
proof  here. 
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From  the  First  Law  of  thermodynamics  in  the  form 


Tds  =  dh  -  vdp 


we  can  immediately  write 


( «•)  =  ( («)  . .  r  (£E, 
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For  isenthalps  within  the  two-phase  region,  however,  ® 

where  ^(p,  T)  =  0  is  the  equation  of  the  vapor -pres sure  curve.  Since 

do  0s  0T 

(•^) .  is  always  positive,  the  slopes  of  the  isenthalps  (‘^).  or  (tt-), 
ax  9  "  V*  u 

must  always  be  negative.  Again,  from  the  First  Law, 
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or 
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The  first  term  on  the  right-hand  side  of  the  above  expression  is  again 
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positive,  as  is  the  third  term,  since  (^).  must  be  negative  in  the 

OX  h 

two-phase  region.  In  the  following  article: 

7.  L.  Riedel,  Kaltetechnik  ^  (1 957) 

it  is  shown  that  the  "logarithmic  curvature"  of  the  vapor -pres  sure 
curve  is  always  negative  for  almost  all  substances  (particularly  for 
hydrogen).  The  exceptions,  such  as  for  water,  occur  only  in  close 
neighborhoods  of  th<.  critical  point;  even  here,  the  validity  of  the 
exception  is  far  from  certain  due  to  the  uncertainty  in  the  critical 
temperature  (this  temperature  can  be  lowered  approximately  one-half 
degree  by  solution  of  quartz  in  water).  Therefore,  the  curvature 

components  (^“l),  and  (  ■^),  =  -  (I”). 

positive. 


must  both  be 


Because  of  such  gross  irregularities,  the  NACA  diagram  (4) 
should  not  be  used  for  design  calculations. 
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The  final  diagram  (5),  an  unpublished  enthalpy-pressure  chart 
due  to  H.  L.  Johnston,  is  believed  by  us  to  be  the  result  of  applying 
an  approximate  method  of  correction  to  Johnston's  data  for  normal 
hydrogen  to  obtain  thermal  properties  of  parahydrogen.  This  approx¬ 
imate  method  will  be  discussed  briefly  below.  The  diagr2Lm,  itself, 
appears  to  satisfy  all  the  conditions  given  above.  In  particular,  the 
terminal  point  of  the  inversion  curve  appears  to  be  only  slightly  above 
27* K.  This  inversion  curve  and  some  of  the  others  referred  to 
previously  are  shown  in  Figure  13. 

It  is  sometimes  claimed  that,  in  the  method  of  calculating 
thermodynamic  functions  from  the  radial  distribution  function, 
statistical  mechanics  obtains  the  result  that  the  derivatives  of  the 
density  with  respect  to  pressure  are  independent  of  the  internal 
molecular  states.  This  "result'*  is  false  in  that  it  is  actually  equiva¬ 
lent  to  an  assumption  introduced  at  the  beginning  of  all  statistical 
mechanical  treatments:  that  the  Hamiltonian  is  separable  into  an 
internal  and  an  external  part.  Nevertheless,  as  an  ad  hoc  assumption, 
relatively  good  approximations  may  be  obtained  by  this  method,  as 
evidenced  by  the  fact  that  at  20*  K  the  densities  of  normal  and  para¬ 
hydrogen  differ  by  only  1/2  percent,  according  to  data  from 

8.  R.  B.  Scott,  Cryogenic  Engineering,  D.  van  Nostrand 

(1959) 

If  h(p,  T)  is  the  enthalpy  at  given  pressure  and  temperature  and 
h*(T)  is  the  enthalpy  of  the  ideal  gas,  the  assumption  of  separability 
of  internal  and  external  states  implies  that  the  quantity  h(p,  T)  -  h*(T) 
is  independent  of  ortho-para  composition.  By  using  the  perfect  gas 
enthalpy  differences  tabulated  in  (2),  approximate  calculations  for 
parahydrogen  can  be  made  from  normal  hydrogen  data:  explicitly, 
the  enthalpy  of  parahydrogen  will  be  given  by  the  approximation 

Vra)''=-^>'  Vormal)'P- 

where  H*(T)  is  the  "heat  of  conversion"  of  the  perfect  gas  from  the 
normal  mixture  to  parahydrogen.  Other  thermal  properties,  such  as 
the  specific  heat  cp  of  parahydrogen  can  also  be  obtained  from  this 
approximation  by  well-known  relations. 

In  conclusion,  use  of  either  the  diagram  (2)  together  with  the 
above  correction,  or  the  unpublished  chart  (5)  (if  available)  should 
give  the  most  accurate  results  in  calculations  with  parahydrogen. 
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PRESSURE  ,  atmospheres 
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However,  either  course  should  be  considered  a  stop- gap  until  the  PVT 
data  presently  being  obtained  directly  on  parahydrogen  at  NBS-CEL  can 
be  used  to  prepare  a  new  temperature-entropy  diagram. 

3.  7  Conclusions 


It  is  planned  to  make  regtilar  visits  to  Convair  in  the  future  to 
provide  assistance  on  Project  Centaur's  cryogenic  engineering  problems 
in  any  way  possible.  In  addition,  visits  are  planned  to  other  instal¬ 
lations  that  might  be  concerned  with  cryogenic  engineering  problems 
relating  to  Projects  Centaur,  Rover  and  Saturn.  It  is  not  expected 
that  the  problems  arising  from  these  visits  will  be  greatly  different 
from  those  already  encotmtered,  however,  such  visits  should  serve 
to  strengthen  coordination  of  the  cryogenic  engineering  effort  on  all 
of  these  projects. 
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Appendix  I 

The  Thermal  Properties  of  Gaseous  and  Liquid  ^drogen  Belov  60°K. 
A  Preliminary  Analysis  of  the  Existent  Data  -  W.  T.  Ziegler 


Summary 


A  literature  survey  covering  the  period  1940-1956,  prepared  by  Mr.  R. 

H.  Kropschot,  has  been  used  to  determine  tdiatis  known  about  the  thermal 
properties  of  liquid  and  gaseous  hydrogen  below  60“K.  !Ilxe  excellent 
summary  prepared  by  WocUey,  Scott  and  Brickwedde  [l]  was  assumed  to 
represent  the  best  available  information  prior  to  about  1944. 

The  literatvire  search  revealed  that  there  was  available  extensive  new 
data,  most  of  it  the  work  of  H.  L.  Johnston  and  his  co-workers  at  Ohio 
State  University.  These  new  data  permit  the  construction  of  P-V-T  emd 
thermodynamic  property  charts  for  the  liquid  region  from  20®  to 
approximately  the  critical  temperature,  55.2®K,  for  pressure  up  to  100 
atm.  The  limited  amount  of  earlier  data  by  Bartholoml[lO]  permit  the  ex¬ 
tension  of  these  charts  to  l4®K. 

Ccraparisons  have  been  made  of  the  compressibility  factors,  Z,  observed 
by  Johnston  et  ed.  with  those  obtained  from  RP1952  [l]  for  several 
isotherms.  In  general  the  ccxnpressibllity  factors  observed  by  Johnston 
were  lower  than  those  calcxilated  from  RP1952  (Table  15)*  The  maxlmuir 
discrepancy  noted  was  about  2  per  cent,  in  the  range  20  -  60®  for  an 
Amagat  density  up  to  500.  In  the  range  below  20  -  55 “K  the  P-V-T  data 
of  Johnston  et  aLl.  is  much  more  extensive  than  the  Leiden  data  available 
at  the  time  RP1952  was  written. 

The  heat  capacity,  C  ,  for  liquid  hydrogen  has  been  shown  to  be  in  good 
agreement  with  the  e^erlmenteLL  data  of  Gutsche  [24]  near  20®K  but  the 
data  of  Gutsche  appears  to  be  approximately  10  -  15  per  cent  too  large 
near  50®K,  in  agreement  with  the  conclusions  expressed  in  RP1952. 

Values  for  C  vs  p  at  constant  T  and  P  have  been  calculated  using  RP1952 
and  the  behailor  of  emd  in  the  criticeO.  region  have  been  discussed. 

The  thermodynamic  properties  ceLLculated  by  Johnston  et  al.  from  the  P-V-T 
data  have  been  compared  with  RP1952. 

The  present  study  is  by  no  means  a  complete  one.  However,  it  is  felt 
that  certain  ideas  on  which  future  study  should  concentrate  have  been 
determined. 
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Introduction 


The  present  study  is  the  result  of  a  project  to  experimentally  determine 
new  data  and  to  re-examine  existing  data  available  on  the  thermal  and 
transport  properties  of  gaseous,  and  liquid  hydrogen  over  the  range  lU- 
60"  and  up  to  approximately  100  atm  pressure.  The  excellent  survey  and 
compilation  of  the  properties  of  hydrogen  by  Woolley,  Scott,  and 
Brickwedde  [  1]  (hereafter  referred  to  as  RP1952)  constituted  the  point 
of  departure  of  the  re-examination,  Mr.  R.  H.  Kropschot  has  made  a 
literature  survey  for  1940-1955  in  ChemiceLl  Abstracts  and  1940-1956  In 
Physics  Abstracts  covering  the  areas  of  thermal  properties,  viscosity 
and  thermal  conductivity.  For  the  years  preceding  1940  the  bibliography 
of  RPI932  has  been  relied  upon,  although  a  number  of  other  earlier  papers 
have  been  examined  for  original  data. 

The  present  study  is  concerned  only  with  the  thermal  properties  of 
liquid  and  gaseous  hydrogen.  In  all  instances  the  word  hydrogen  Im¬ 
plies  *  normal**  hydrogen  (75  per  cent  ortho,  25  per  cent  para)  unless 
otherwise  stiyl?  nted.  In  the  earlier  work  on  hydrogen  the  ortho-para 
composition  is  vinkncwT»;  it  is  assumed  to  be  for  normal  hydrogen. 

The  properties  which  are  discussed  have  been  grouped  under  the  follow¬ 
ing  broad  headings: 

P-V-T  relations  for  gaseous  and  liquid  hydrogen 

Vapor  pressure 

Critical  constants 

Joule-Thomson  effect 

Heat  capacity 

Velocity  of  sovmd  and  relation  to  C  /C 
Thermodynamic  properties  ^  ^ 

A  more  detailed  listing  is  given  in  the  Table  of  Contents. 

Time  available  has  not  permitted  a  critlced.  analysis  of  the  data  avail¬ 
able  in  these  areas.  However,  an  atten5>t  has  been  made  to  detennine  the 
probable  effect  of  the  newer  data  on  the  conclusions  reached  in  RP1932, 
to  indicate  ^ere  new  calculations  are  likely  to  be  useful,  and  to 
determine  where  new  experimental  studies  should  be  undertaken. 

)bich  of  the  data  on  the  thermal  properties  of  hydrogen  which  has 
appeared  since  RP1932  was  con5)leted  are  the  result  of  the  extensive 
Investigations  of  H.  L.  Johnston  and  his  coworkers  at  Ohio  State 
University.  Much  of  this  work  has  been  reported  in  Technical  Reports 
rather  than  in  the  open  literature.  In  other  reports  little  effort 
has  been  made  by  the  authors  to  anedyze  the  results  critically  with 
respect  to  existing  data.  The  work  of  Johnston,  et  al.  represents 
a  considerable  extension  of  the  data  available  In  '  /42-46  when  RP1952 
was  being  written.  Thus,  P-V-T  data  for  the  liquid  in  the  range  20- 
JJ'K  for  pressures  of  2-100  atmospheres,  P-V-T  data  for  gaseous  hydrogen 
below  33“K,  the  heat  capacity  from  the  boiling  to  the  critical  point 
constitute  information  which  for  the  most  part  was  unknown  in  1942-46. 
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Cbe  of  the  prohlems  which  arises  in  the  combination  of  data  from  various 

sources  is  that  of  the  absolute  temperature  scale  used  by  the  different 

investigators.  The  temperature  scale  used  throughout  the  present  paper 

has  for  its  ice  point  T  =  273*16®  (the  NBS  scale  of  RP1952).  Much  of 

the  Leiden  data  is  given  in  the  absolute  centigrade  scale  with  T  =  275 *09*^. 

V/here  necessary,  the  corresponding  temperatxire  on  the  scale  with 

T  =  273*16®K  was  obtained  from  the  relation, 
o 

^NBS  "  (273.09  +  t®c) 

This  method  of  conversion  was  used  by  Woolley  et  elI  [l]  (conversation 
with  R.  B.  Scott).  All  temperatures  unless  other  specified  refer  to 
®K  on  this  scale. 

P-V-T  Relations 


a.  Gas 

The  P-V-T  relations  for  gaseous  hydrogen  prior  to  about  19^46  have  been 
summarized  in  Fig.  6  of  Reference  [l]  and  used  to  construct  the  Z  vs  P , 

T  table  of  this  work  (RP1952).  Examination  of  Fig.  6  shows  that  all 
these  data  were  the  result  of  measurements  made  at  Leiden.  The  original 
Leiden  data  below  60.55“K  has  been  reexamined  (6  land  a  sumnary  of  it 
collected  in  Table  I. 

Since  1946  a  large  amount  of  additional  data  in  the  fom  of  approximately 
50  isothermb  extending  from  20  to  300®K  and  1  to  200  atmospheres  has 
appeared  [  7^  10  ]  .  Fifteen  of  these  isotherms  are  for  tenperatures 

in  the  range  20-60®K,  seven  of  them  being  below  the  critical  tenperature 
03.24®K),  The  data  are  stunmarized  in  Table  II.  This  new  work  is  the 
result  of  investigations  of  H.  L,  Johnston  and  coworkers  at  Ohio  State 
University.  It  will  be  referred  to  hereafter  as  by  Johnston.  Johnston 
[8  lhas  fitted  vlrial  equations  of  the  form 

g  =  z  =  i+  |  +  -I2+  ^ 

to  his  isotherms  for  T  ^  55.10®K  and  compared  Z  ”  ^  oi  over  the  range 
of  his  data.  He  states  [8,  p  4  Ithat "  it  is  ev?dent  t£§t^the  representation 
of  these  data  by  the  vlrial  equation  yields  results  which,  on  the  average, 
are  In  agreement  with  the  e]q>erlmentsl  data  to  one  part  In  1000,  or 
better,  except  at  the  lowest  temperatures.®  At  53.96®K  It  was  found 
Impossible  to  fit  the  data  to  the  six-term  vlrial  equatlm.  For 
temperature  T^  33.00®K  Johnston  [7  Juses  fewer  vlrial  coefficients  and 
does  not  explicitly  give  the  comparison  ^  ^  -  Z  _  .  In  none  of  these 

papers  is  a  coanparlson  made  between  the  Leiden  dam,  the  Z  vs  P  ,  T  data 
In  RPI932  and  the  Johnston  results.  However,  it  is  clear  that  the 
Johnston  data  is  much  more  extensive  in  the  razxge  T<  ~3^®K  than  is  the 
Leiden  data. 

A  few  Isotherms  very  near  the  critical  point  have  also  been  measured  by 
Hoge  and  Lassiter  [8  1 . 
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TABJ.E  T 


Isotherm 

Range  of  Amagat 

Points 

Leiden  Conm 

t°C 

T°K 

density, 

In  data 

No. 

(RPI952) 

-212.82 

60.55 

74-194 

6 

100a 

157-525 

5 

9Ta,99a 

70-90 

3 

165b 

-217 

55.77 

78-195 

6 

100a 

285-559 

3 

97a,  99a 

70-90 

7 

165b 

202-544 

7 

172a 

-225.37 

47.74 

69-90 

5 

165b 

224-455 

8 

172a 

16-25 

4 

l88e 

-251.38 

41.70 

72-89 

3 

165b 

246-533 

l4 

172a 

16-25 

4 

l88e 

-236.56 

16-25 

5 

l88e 

-236.29 

36.79 

75-89 

4 

165b 

268-602 

15 

172a 

-238.29 

54.81 

80-88 

5 

165b 

285- 60 

12 

172a 

-259.90 

35.19 

79-89 

5 

165b 

298-700 

22 

172a 

-24i .87 

51.22 

79-89 

5 

165b 

16-24 

4 

l88e 

-245.88 

29.21 

80-90 

7 

165b 

-248.32 

24.77 

19-25 

4 

l88e 

38 


TABLE  II 


Isothem 

T  ’K 

Pressure  range 
atm 

Aoagat  density^ 
p«  range 

Reference 

20.59 

0.20  -  0.99 

2.8  -  14.5 

I7I 

22.58 

o.to  -  1.67 

5.1  -  25.2 

7 

2U.65 

0.26  -  2.6 

3.0  -  35.1 

7 

26.T5 

0.51  -  U.o 

5.3  -  53.1 

7 

28.85 

0.58  -  6.0 

5.6  -  82.5 

7 

50.86 

0.1».9  -  8. *1^5 

4.4  -  120.8 

7 

55.00 

55.2MT^) 

0.98  -  11.5 

8.4  -  178.0 

7 

33.96 

1^.0  -  1^.2 

44.9  -  859 

8 

35.10 

0.98  -  14.5 

7.8  -  210 

7 

16.2  -  190.7 

567.2  -  843 

8 

37.61 

5.4  -  199.5 

6.6  -  190 

8 

4o.09 

0.95  -  17.9 

6.6  -  190 

7 

23.8  -  200.1 

300.5  -  820 

8 

U5.10 

0.9  -  17.2 

8.1  -  130.6 

7 

22.8  -  202.8 

190  -  783 

8 

50.09 

2.06  -  16 -.2 

11.4  -  102.3 

7 

22.8  -  200 

153  -  7^*4 

8 

55.09 

0.97  -  16.6 

4.8  -  91.7 

7 

25.6  -  198.6 

136  -  704 

8 

60.05 

0.93  -  30.5 

4.2  -  158 

7 

59.8  -  195.7 

211-665 

8 

65.96 

17.5  -  205 

80.7  -  647 

8 
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Interccjanparlson  of  Johnston  data  with  RP1932  and  Leiden  data.  Ihls 
Interccaiqparlson  was  carried  out  at  selected  temperatures.  Mid  presstires 
in  order  to  detennlne  how  closely  the  new  Johnston  data  fitted  the 
Z  vs  p,  T  table  given  in  RP1932.  Certain  isotherms  taken  frcm  the 
origined  Leiden  data  were  also  compared  with  the  Table  15  in  RP1932 
to  detennlne  how  well  this  Table  reproduced  the  orlgincd  Leiden  data 
on  which  it  was  based.  All  comparisons  have  been  made  below  *'60 "K. 

In  order  to  make  this  comparison  the  Amagat  density  corresponding 
to  ftll  points  for  T  ^63.96“  of  the  Johnston  data  have  been  calculated 
using  the  relation 


P 

Z 


P  V 

TWJ 

PV 

RT 


(2) 

(3) 


where  P  and  PV  are  given  by  HLJ  in  atm  and  cc.  atm  mole  ”  ,  respectively. 
The  value  of  R  used  by  Johnston  was  R  =  82.05667  cm5  atm  deg“l  [8,  p.  3], 
The  value  of  Vq  [the  molal  volume  of  0“C  and  1  atm  pressure]  used 
was  22U28.5  cm^  mole  taken  from  RP1932,  p.  396.  Z  has  been 
calculated  by  Eq.  (3)  where  this  had  not  been  done  by  Johnston. 


The  Intercomparisons  have  been  carried  out  in  two  ways.  In  the  hlgb 
density  range,  p  >  2CX),  T  <  56“K,  the  comparison  was  made  by  four-point 
Lagranglan  interpolation  in  the  Table  13  of  RP1932.  In  the  low  density 
region,  p  <  200,  T  <  56 “K,  the  Intercon^arlson  has  been  carried  out  by 
means  of  the  relation 

(4) 

where  A  and  C  have  been  taken  from  Table  19  of  RP1932.  Eq.  (!(-)  provided 
a  simple  means  of  obtaining  FV  since  the  Table  13  1b  calculated  from 

RT 

Eq.  (1|-)  together  with  Table  I9  [see  remarks  in  RP1932,  p.  ^30]. 


High  density  range,  p  >  200.  Z  has  been  calculated  from  Table  13 
for  each  point  on  the  Johnston  Isotherms  T  =  l4-5*10‘’,  37.61  ®K,  as  a 
function  of  p  and  the  difference  -  Z  .  ^  plotted  in  Fig.  1.  A 

similar  comparison  has  been  made  for  Leiden® Isotherms  T  =  47.7^"  and 
T  =  36.79“  lii  order  to  examine  the  agreement  between  Table  I3  and  the 
Leiden  data  used  to  construct  the  Table  13  in  this  range.  [Data  for 
the  43.10”K  Isotheim  of  Johnston  in  the  range  p  <  200  Is  Included.  It 
was  calculated  using  equation  (4)  as  described  below] . 

Examination  of  Fig.  1  shows  that  Z__g  -  Z  .  ^  is  less  than  2  x  lO"^  for 
all  points  except  one  of  the  Lelden^^ta  ^ereeis  the  Jcdinston  data 
gives  Z_^  -  Z  ,  g  as  large  as  6  x  10"5.  There  seems  to  be  more 
scatter”u  the^Johnston  data,  but  also  some  indication  of  a  trend  for 
(HLJ)  to  be  systematically  smaller  at  hl^  p  than  the  Leiden  data. 
Instance  at  T  =  57.6l*K  p  =  474.16  =  0.4152  Z  .  (HLJ)  = 

0.40920,  the  Z  observed  by  HLJ  being  1.45”^r  cent  low.  °  ® 


Figure  1 

Comparison  of  Z^-Z^iden 

approx,  the  same  isothems 
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Lew  density  range,  p  <  200  ,  T<^"K.  been  calculated 

as  a  function  of  p  for  all  points  given  on  the  JMnston  Isotherms  at 
T  =  20.59“,  24.65%  26.75“,  28.85“,  50.86“,  55.00"  (all  below  T  =  55.24") 
and  55.10"  and  45.10"  using  equation  (4).  Appropriate  values  o?  A  and  C 
for  these  temperatures  were  obtained  by  plotting  the  values  of  A  and  C 
given  In  Table  19  (RP1952).  The  difference  -  Z  .  ^  has  been  plotted 

vs  p  In  Figure  2  for  all  Isotherms.  A  similar  coo^wr  is  on  Is  shown  for 
the  Leiden  Isotherms  at  T  =  29.21",  51.22",  55.19“,  and  54.81 "K  In  Figure  5* 

It  will  be  noted  from  Figure  5  that  Z_^  “  ^  quite  small  for  T  =  54.81" 

and  55.19“.  The  agreement  is  less  satisfactory  for  T  =  31.22"  and  29.21"K. 
This  behavior  Is  In  agreement  with  conclusions  to  be  drawn  from  Figure  9 
(RPI952)  p.  429,  where  the  values  of  A  and  for  the  first  two  isotherms 
fall  on  the  curve  used  to  obtain  Table  I9,  while  the  experimental  values 
of  A  for  T  =  51.22  and  29.21"  fall  off  the  smooth  curve  drawn,  hence  one 
would  expect  some  disagreement  between  Table  I5  and  the  Leiden  data  in 
this  temperature  range .  It  should  be  noted  that  these  four  Isotherms 
represent  a  very  limited  range  of  p,  whereas  the  Johnston  data  cover  a 
much  wider  range  of  p.  Examination  of  Figures  6  and  8  of  RPI952  reveals 
how  sparse  the  Leiden  data  really  are  in  the  temperattire  range  20-56" 
for  p  <  200 . 


Comparison  of  the  Johnston  data  for  Z  with  that  calculated  from  RPI952, 
as  shown  in  Figure  2,  shows  that  as  one  goes  from  20.59“  to  45.10"  the 
agreement  between  the  two  is  good  at  T  =  20.59“  where  p  Is  riiih.11  (and 
limited  by  the  saturation  pressure),  gets  progressively  poorer  as  T 
Increases  to  about  55 “K  and  then  Improves  again.  We  may  note  that  Z 
(for  a  given  condition)  of  Johnston  is  almost  invariably  smaller  than 
the  corresponding  Leiden  value. 


The  percentage  discrepancy  between  the  Johnston  data  and  the  NBS  Table 
13  of  RPI952  varies  with  the  value  of  Z.  At  T  =  35.00"K  where  Z  «  0.5 
Z_^  is  1.85  per  cent  smaller  than  for  p  =  I78.  This  is  the  maximum 
dil^repancy  noted  except  for  isolatea  points . 


ISie  possible  effect  of  an  error  in  the  temperature  measurement  of  the  Isotherm 
was  examined  at  35.00“K.  Decreasing  the  temperature  to  32.90“  Improved  the 
fit  somewhat.  It  appears  that  the  teii5>erature  would  have  to  be  decreased 
to  about  32.6"  to  give  the  values  of  Zcalc  from  the  NBS  Table.  Such  an 
error  in  the  temperat\ire  measurement  seems  highly  unlikely  in  this 
ten^ratxire  range. 

One  cooiparlson  of  Johnston  s  data  with  his  vlrlal  equation  for  T  =  30.86" 
and  P<  ^saturation  n«de.  The  equation  used  as  [7]  . 

PV  1  +  B+  C_  +  D 
RT  -  V  ^  v5 

where  B  =  -  80.75  cc/mole 

C  =  2078  (cc/mole)^ 

D  =  50,000  (cc/mole)^ 

This  equation  was  found  to  fit  the  experimental  data  to  better  than  2  parts 
per  1000  except  at  one  point  in  the  range  p  *  120.8.  This  fit  is  as  good  as 
can  be  obtained  with  equation  (4) . 

It  seems  likely  tliat  the  fit  for  other  temperatures  below  33 *K  will  be  equally 

good. 
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b.  Liquid 

The  P-V-T  relations  for  liquid  hydrogen  (n  -  H„)  In  Reference  1  (RP1952) 
are  based  on  the  experimental  measurements  of  Mathias,  Crcmmelln,  and 
Kamerlln^  Onnes  Leiden  Comm.  No.  1546b  ll2  ]  for  the  saturated  liquid 
over  the  range  20  -  53 and  upon  the  measurements  of  Bartholcni^  110  1 
for  three  Isotherms  having  the  following  range  of  data: 


Pressure  range 
range _ 


Freezing  pressure* 
atm 


16.45  8  -  82  81.7 
18.24  5  -  157.5  148.5 
20.35  1  -  236.5  232.9 


*Calc.  from  Equation  7.l8,RP1932. 

Since  the  appearance  of  RP1952>  Johnston  and  his  coworkers  [9,  35  1  have 
measured  the  P-V-T  relations  for  9  Isotherms  covering  the  range  20.34* 
to  32.58*K  from  approximately  10  to  100  atmospheres,  niese  results 
have  been  presented  in  the  form  of  P  vs  T  plots  at  constant  V,cuid  P  vs  V 
plottat  consteint  T[  9]  •  No  new  data  were  obtained  for  the  volume 
of  the  saturated  liquid;  Johnston  et  al  extrapolated  their  data  to  the 
satiiration  volirmes  given  by  Mathias  et  al  cited  above.  Ccxnparlsons  of 
Johnston's  data  for  T  -  20.38“  with  that  of  Ba3rtholcm6  at  T  =  20,59“ 
shows  that  the  agreement  is  1  per  cent  or  better  over  the  range  9  * 

106  atm  (where  ccn^jarison  is  possible),  Johnston  et  al  [351  have 
extrapolated  their  volxime  data  to  I50  atm  and  have  used  their  smoothed 
P-V-T  data  to  calculate  the  thermodynamic  properties  of  the  liquid 
relative  to  the  properties  of  the  saturated  liquid. 

Friedman  and  Hllsenrath[  5]  have  used  Johnston's  data  [  9  1  to  prepare 
a  table  of  T  vs  P  values  at  constant  V  and  compressibility  factor  Z 
vs  P  at  constant  V  over  the  range  I5  -  55‘’K.  In  a  subsequent  publlcatl(xx 
Friedman  [  4  ]  has  prepared  a  table  of  Z  vs  T  at  constant  P  over  the 
range  I6  -  35“K.  In  both  instances  the  pressure  goes  to  120  atm. 


c .  Saturated  vapor  and  liquid 


No  new  data  for  the  saturated  vapor  or  liquid  has  appeared  since  the 
summary  presented  in  RP  1932.  It  appears  that  the  data  given  for 
n  -  Hp  (as  liquid)  in  Table  51  l460of  RP1952  should  be  adequate.  It 
should  be  noted  that  the  data  given  for  p  -  in  Table  51  differ 
slightly  from  those  for  n  - 

Vapor  Pressure  Relations 


White,  Friedman,  and  Johnston  [  15I  have  recently  measured  the  vapor 
pressure  of  n  -  from  the  boiling  point  to  the  critical  point.  In 

another  paper  they  have  measured  T  and  P  by  the  visual  observation 
of  the  disappearance  of  the  meniscus  and  found  T^  _  33 .24^  ®K  and 
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P  =  12.797  atm  [I9]  • 
c 

They  have  expressed  their  measurements  In  the  form  of  an  equation 

(5) 

Which  gives  20.U0“K  for  the  normal  b.p.  and  P  =  12.77*  atm  for  T  =  53.2Uj^“K. 
Comparison  of  this  equation  was  made  with  one  given  In  RP1932.  %t 
Identical  pressures  the  temperature  calculated  from  each  equatlcm  was 
compared.  The  difference  In  temperatxire  varied  approximately 

linearly  with  T  reaching  +0.18“  at  the  critical  t^perature.  It  should 
be  pointed  out  that  the  equation  from  RP1932  Is  actually  based  on  the 
range  lU  -  20*K,  but  was  stated  (RP1932  p  ^5^)  to  represent  the  Leiden 
data  to  the  critical  temperature  (33.19*K)  within  the  limits  of 
experimental  accuracy. 

A  detailed  comparison  of  the  Johnston  and  Leiden  data  should  be  made  to 
determine  the  nature  of  the  agreement  directly.  GrlUy  [l6]  heus  also 
reported  the  vapor  presstire  of  n  -  Hg  to  3  atm  pressure. 

Hoge  and  Arnold  [l^l-]  have  recently  measured  the  vapor  pressures  of 
p  -  Hg  (99*71  per  cent  p,0.21  per  cent  ortho),  Dg  and  HD  from  the 
normal  b.p.  to  the  critical  point.  Their  experimental  results  have 
been  e^qpressed  for  p  -  Hg  in  the  foxm  of  a  table  permitting  linear 
Interpolation  except  above  30 “K.  It  Is  Interesting  to  note  that  Hoge 
and  Lassiter  [l8]  found  the  critical  constants  of  p  -  Hp  to  be  T  » 

32.994%  P  =12.770  atm,  V  =65.5  cc/mole.  ^ 

It  should  be  noted  that  the  vapor  pressure  of  n  -  Hg  Is  sll^tly  less 
at  the  same  tenperature  than  that  of  p  -  Hg^Csee  BPI932)  for  the  range 
l4  -  20®K.  Johnston's  equation &r  n  -  Hg  fEquatlon  5  above]  has 
this  property  at  20.4*  and  32.99^"K.  In  should  be  noted  that  o-p  Hg 
solutions  do  not  follow  the  Ideal  solution  laws  [RP1932,  Reference 

14,  15,  17]  . 

Critical  Constants 
n  -  ^ 


Woolley  et  al  (RP1932)  have  suonarlzed  the  critical  constant  data  for 
n-Hg.  The  latest  Leiden  values  as  reported  by  Woolley  et  al  are 

T  =  33.19“  (on  the  basis  T  =  273*16“)  P  =  12.751  atm  and  =  3**^. 
(&nsgat  density).  However, °on  the  basis  of  the  P-V-T  data  Woolley  et  al 
finally  arrive  at  a  compromise  choice  as  follows: 

T  =  33.19“  P  =  12.98  atm  p  =335  V  =  66.95  cc/mole  =  0.3191 

c  c  c  c  c 

More  recently  Johnston  et  al  [19I  have  made  a  direct  determination  of  the 

critical  point  by  observing  the  disappearance  of  the  meniscus.  Hiey  report 

T  =  33*24.  *K  P  =  12.797  atm. 
c  4  c 
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p  -  H2 

Recently  Hoge  and  Lassiter  [l8]  have  measured  the  critical  constants 
of  p  -  Hg  (99.71  per  cent  p,  0.21  per  cent  o),  HD  and  e-D2 
Isothenn^method .  Their  results  for  p  -  Hg  are 

=  32.99»v"K  =  12.770  atm  =  65.5  cc/mole  =  O.309 

0^  = 

Joule -Thomson  Effect 

Johnston  et  al  [20]  have  measured  the  Joule-Thcmson  effect  in  normal 
Hp  directly  from  and  at  room  tenq^rature .  Most  of  the  measurements 

vere  made  In  the  range  -  dO‘’K.  The  results  of  these  measurements 
have  been  used  by  them  [21]  to  calculate  the  thermodynamic  quantities 
H-H®,  C  -  C  “  and  S-S“  over  the  range  65“  -  300®  and  at  even  pressures 
from  5  “  200^atm. 

Baehr  has  used  the  same  data  to  calculate  the  quantities 


and  Cp  over  the  range  60®  to  85®  from  P  =  0  to  P  =  I80. 

Both  of  these  latter  papers  will  be  discussed  In  connection  with  the 
heat  capacity  and  thermodynamic  property  section  of  this  repoirt. 


Heat  Cat 


The  heat  capacity  of  a  substance  may  be  expressed  In  a  number  of  we^s 
depending  upon  the  partlcvilar  process  Involved  or  the  method  of  calcu¬ 
lation  used  to  obtain  it.  In  this  section  we  will  discuss  the  follow¬ 
ing  both  for  gaseous  and  liquid  hydrogen: 


C  C  r  X  C  _  C  /C  ,  C  -  C  ". 
p  '  p(satd.},  satd.,  p'  v'  v  v 


Woolley  et  al  (RP1932)  have  svonmarlzed  the  e;q)erlmental  data  which 
consists  of  C  meastirements  over  the  range  I6.5  -  37*K  and  at  seven 
pressures  (ll^-  99  kg/cm^)  by  Gutsche  [2k]  and  C  measurements  by 
Eucken  [26]  and  Eartholome  and  Eucken  [23]  In  the  solid  and  liquid 
states.  No  new  measurements  have  appeared  for  normal  hydrogen. 
Johnston  et  al  [23]  have  measured  C  ^  for  liquid  parahydrogen  from 
12.7®  -  20.3®K  and  In  a  later  paper^S^e  reported  for  liquid 

p-Hg  from  the  boiling  point  to  the  critical  point  [22], 


Woolley  et  al  (RP1932)  have  used  their  equation  of  state  table,  together 
with  the  experimental  C  data  of  Gutsche,  to  construct  a  plot  of  C  vs  T 
for  several  constant  pi^ssTires  for  coopressed  gaseous  and  liquid  hydrogen 
[their  Fig.  27].  This  Figure  covers  the  range  15  -  42 ®K.  Chelton  et 
al[5]have  used  the  equation  of  state  data  of  RP1932  to  calculate  C 
for  gaseous  hydrogen  over  the  range  20  -  300®K  at  constant  pressures 
of  12.98  (P  )  15,  20,  30,  40,  60,  100  atms.  (Also,  1,  4.84-  and  10  atm.) 
The  data  of^Gutsche  have  been  used  to  conq>lete  the  Chart  In  l^e  highly 
compressed  fluid  (liquid)  region. 
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Uie  extensive  P-V-T  measurements  of  Johnston  and  coworkers,  together 
with  the  Joule-Thorns on  measurements  of  this  group,  and  their  recent 
heat  capacity  measvirements  on  parahydrogen  (liquid)  permit  one  to 
obtain  C  by  making  use  of  appropriate  thermo-dynamic  relations.  On 
the  one  Sand,  these  calcvilations  permit  an  independent  check  on  the 
experimental  measurements  of  Gutsche  and,  on  the  other,  an  independent 
check  or  conparison  with  the  C  results  of  Chelton,  et  al.  [^1  for  the 
gas.  An  analysis  of  the  curves  of  Chelton,  et  al.  will  be  made  from 
this  point  of  view  and  the  data  of  Gutsche  ccmpcured  with  the  calculations. 

Some  connents  will  also  be  made  on  the  direct  measurement  of  C  and 
C  ,  the  behavior  of  C  and  C  near  the  critical  region,  the  meSsurement 
olF  C  and  C  from  velBclty  o¥  sound  measurements  and  the  direct  meas\ire- 
ment^of  the^C  /C  ratio. 

p'  V 

a.  of  gaseous  n-Hg 


1.  Calculation  from  Jcnle-ThonuBon  measurements 

Vfhlte  and  Johnston  [  3^  have  \ised  the  earlier  Joule-Thomson  measurements 
of  White* Friedman  and  Johnston  [  20]  from  63°  -  ^OO^K  and  1  -  200  atm 
to  calculate  H  -  H®  at  5“  Intervals  from  65“  to  300“K  and  for  17 
constant  pressures  from  2.5  to  200  atm.  Shis  Thble  was  used  to  obtain 
C  -  C  ®  by  linear  interpolation.  Combining  these  values  of  C  -  C  * 
wSth  tSe  values  for  C  ®  given  by  Hllsenrath  et  al  [  2  ]  one  can  ^  ^ 

obtain  C  and  coopare^the  results  with  the  calculation  of  C  made  by 
Chelton  it  al.  ^ 


Baehr  [  2l]  has  6Q.S0  used  the  Johnston  J.  T.  data  I  20]  to  compute  C 
over  the  range  60  -  85 "K  and  o  to  I80  atm. 
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Conparison  of  the  results  of  Johnston  for  C  in  the  range  65®  -  100“K 
was  made  at  40,  60,  and  100  atm  with  the  corresponding  C  results  taken 
from  the  Chart  of  Chelton  et  al  [  5  ]  .  Down  to  75“K  the  ¥esultB  agree 
to  within  about  3  per  cent  (0.2  cal/mole  deg).  At  TO®  the  difference 
is  7*3  per  cent  for  100  atm,  3*8  per  cent  at  60  atm  and  2  per  cent  at 
Uo  atm,  all  lower  than  the  Chelton  values.  The  C  values  obtained  by 
Baehr  agree  less  well  with  the  Chelton  results  tl£i  do  those  of  Johnston. 


It  appears,  therefore,  that  these  J-T  results  of  Johnston  et  al  agree 
with  the  Chelton  results  reasonably  well  to  about  70®K. 

2.  CalcrOations  from  Equation  of  State  Data 

Johnston  and  coworkers  [3^]  have  used  their  extensive  P-V-T  measure¬ 
ments  for  gaseous  normal  hydrogen  to  obtain  C  -  C  ®  over  the  range 
20  -  300®K  for  pressures  from  o  to  200  atm.  ftie  rlsuJ.ts  are  presented 
in  the  form  of  a  graph  of  C  -  C  ®  vs  P  for  a  series  of  isotherms.  In 
the  region  of  Interest  of  toIs  riport  the  Isotherms  are  26,  28,  30,  32, 
36,  38*  40,  45,  50,  55^  60,  65  and  at  5®  intervals  to  100®K  (there¬ 
after  at  10®  Intervals  to  300®).  Unforttinately  no  tabular  data  is 
presented  so  that  C  -  Cp®  must  be  read  from  a  graph. 
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These  calculations  were  made  by  machine  computation  using  the  vlrlal 
eqviatlon  data  to  represent  the  P-V-T  relations. 

In  the  range  70  -  100°K  and  for  pressures  of  4o,  60,  and  100  atm 
the  values  of  C  calc\ilated  from  the  C  -  C  ®  graph  of  Johnston  agree 
within  about  3  fer  cent  or  better  with^the  ?alue  derived  by  Chelton  et 
al  from  the  NBS  P-V-T  relations  In  RP193S*  In  the  range  50  -  70“  the 
agreement  is  within  6  per  cent.  However,  the  error  In  reading  the 
Johnston  plot[  Fig.  2  of  Reference  3^1  nay  be  as  much  as  3  per  cent. 

For  temperatures  less  than  50“K  the  agreement  becomes  progressively 
poorer.  In  fact,  for  T  <  —42“,  C  -  C  “  becomes  negative  for  pressures 
greater  than  40-60  atmospheres.  ^  ^  This  is  a  clear  indication  of 

error  In  C  -  C  “.  In  view  of  the  general  agreement  observed  between 
the  JohnstBn  anS  Leiden  P-V-T  data  and  the  NBS  P-V-T  table  In  this 
region  It  seems  probable  that  the  reason  for  the  discrepsmcy  lies  not 
In  the  P-V-T  data  itself,  but  In  the  failure  of  the  vlrlal  equations  to 
properly  fit  the  P-V-T  Isotherm  data.  White  and  Johnston  make  no 
specific  comments  about  this  pecularlty  in  C  -  C  “  but  do  remark  that 
“the  calc\ilated  thermodynamic  properties  exh?biteB  peculiarities  arising 
from  vincertalntles  in  the  derivatives  of  the  vlrlal  coefficients  [  34-  ]  .* 
In  another  place  It  is  remarked  by  Johnston  [  8  ]  that  the  representation 
of  the  compressibility  factor  data  by  the  vlrlal  equations  yields  results 
which,  on  the  average,  are  In  agreement  with  the  experimental  data  to 
one  part  in  1000  or  better,  except  at  the  lowest  tenperatvures.  “At 
33«96“K  It  was  Impossible  to  fit  the  data  with  a  six- term  vlrlal 
equation.  L8]  * 

For  temperatures  and  pressures  below  T  ar  35.24“K  and  P<  P  =12.8  atm 
we  have  earlier  pointed  out  feee  Section  on  P-V-T  relations^  }that  the 
Johnston  P-V-T  data  is  much  more  extensive  than  the  Leiden  data  on  which 
the  Z  vs  p,  T  chart  of  RP1932  is  based.  Furthermore,  there  is  a  small 
but  cleeur  difference  between  the  values  of  Z  found  by  Johnston  and  those 
read  from  Table  13  of  NBS  RPI932  for  the  same  p  and  T .  This  difference 
In  P-V-T  data  may  be  expected  to  give  C  -  C  “  values  which  differ 
among  themselves  depending  on  the  data  Bsed.^  Ihe  C  vs  T  charts  given 
by  RPI932  [  page  4o4  ]  and  by  Chelton  et  al  [  5  ^  show^lines  only  for 

C  vapor  and  C  for  p  =  5  kg/cm^  (4-. 84  atm),  both  calculated  from 
psatd  P 

the  P-V-T  data  in  ^1932. 


b. 


for  liquid  n-Hg 


Heat  capacity  data  for  n-Hg  in  the  liquid  state  are  rather  sparse.  The 
heat  capacity  of  the  saturated  liquid,  C  ,  has  been  summarized  for  the 
range  l4  -  20“K  by  Woolley  et  al  [  1  ] in  Iheir  Table  37.  No  measurements 
of  C  for  T  >  — 22“K  appear  to  have  been  made.  Johnston  and  his  coworkers 
have^measured  C  for  liquid  p-Hp  over  the  range  l4  -  20.3“K  [23]  and 
frcm  18.28“  to  3i.49“K  [22]  .  The  agreement  between  the  values  of  C 
found  by  Johnston  for  p-Hp  and  those  given  in  RP1932  for  n-Hp  is  very 
good  over  the  range  l4  -  20“K.  In  the  absence  of  data  for  n^g  above 
20“K  and  because  of  the  great  similarity  to  be  expected  between  n-Hg 
and  p-Hp  [the  critical  constants  eind  vapor  pressure  curves  are  very 
nearly  the  same]  it  will  be  assiuned  that  C  for  n-Hp  is  the  same  at 

each  temperature  as  that  found  by  Johnston  et  al  for  p-Hp  ovarthe 
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range  20  -  52 °K. 

Johnston  et  al  [  55  1  have  used  their  earlier  equation  of  state  data  for 
liquid  n-H^  [  9  1  to  calculate  a  table  of  values  of  C  -  C  for  P  =  2 

atm,  to  P  =  150  atm  at  one  degree  Intervals  from  ^  ^satd  t  =  20®K 
to  T  =  55'’K.  In  order  to  use  this  table  to  compute  C  it  Is  necessary 

to  have  both  C  and  C  .  Now  the  relation  between^these  quantities 

Psatd  ^ 

is  given  by  the  relation  [  see  for  example  Dodge,  Chem.  Eng.  Themodynamics 
p.  252  ff  ] 

C  -  C  =  -  T  ^ 

®  Psatd  (.jtJ  p  g  (6) 


where  evaluated  for  the  liquid  at  T  =  p  =  and  ^  = 

slope  of  the  vapor  pressure  equation  at  T  ,  P  .  Johnston,  et  al.,  [35] 
also  give  a  table  of  the  quantity 


Thus,  by  canblnlng  the  tabulated  values  given  by  Johnston,  et  al  [35] 

1  /dV\  dP 

for  V,  ^  \  5^ y  f  ^  calculated  from  their  vapor 

P  ^  ^satd 

pressure  equation  for  n-Hp  [15 1  one  can  calculate  for  liquid  hydrogen 

as  a  function  of  P  and  T.  The  calculation  is  sizramarlzed  in  the  following 
Table  III. 
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TABLE  III 


T„^ 

P(8atd) 

atm 

dP 

df 

atm 

aeg 

V(8atd) 

cc/mole 

l/dV\ 

V  V5f/p 

deg-l 

C 

8 

c  -  c 

®  ^satd 

Cal/mole  'K 

C 

^satd 

20 

0.90908 

0.2798 

28.27 

0.0159 

4.50 

-0.218 

4.72 

22 

1.5583 

0.4035 

29.24 

0.0186 

5.11 

-0.290 

5.4o 

24 

2.5178 

0.5613 

30.59 

0.0233 

5.68 

-0.412 

6.09 

26 

3.8243 

0.7512 

31.94 

0.0311 

6.45 

-0.625 

7.08 

28 

5.5480 

0.9801 

33.91 

0.0430 

7.TO 

-0.989 

8.69 

30 

7.7775 

1.2593 

37.03 

0.0720 

9.80 

-1.957 

11.74 

32 

10.628 

l.6o44 

18 

53 

35.24 


Cbe  may  conqiare  the  experimental  results  of  Gutsche  [24] 
obtained  frcm  the  C  -  C  tables  of  Johnston,  et  al. 

P  Psatd 

results  are  con5>ared  In  the  following  Table  IV. 


with  those 
At  20“K  the 
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TABLE  IV 


T  =  20*K 


p 

kg/cm^ 

P 

atm 

c  -  c 

®  ®satd 

C 

P 

cal /mole  *K 

C 

P 

Calc,  from 

HLJ  [35 1 

Gutsche 
(smooth  curve 

0.909  (V,P) 

0 

4.72 

11 

10.65 

-0.17 

4.55 

4.60 

24 

23.23 

-0.36 

4.36 

4.46 

40 

38.72 

-0.54 

4.18 

4.15 

59 

57.1 

-0.72 

4.00 

3.96 

79 

76.46 

-0.87 

3.85 

3.84 

99 

95.8 

-1.00 

5.72 

3.68 

These  results  are  in  remrkably  good  agreement,  considering  the  fact  that 
Gutsche  estimates  his  results  as  accurate  to  ±  3  per  cent.  Slmileur 
comparisons  have  been  made  hy  me  at  24*,  26*,  28*  and  ^0*K.  The  agree* 
ment  Is  poorer  at  the  higher  temperatures,  the  C  values  calculated  from 
the  Johnston  Table  being  lover  than  measixred  by  Gutsche  by  as  much  as 
1^  per  cent.  The  results  at  T  s  30*K  axe  given  In  the  foUovlng  Table 
V. 


TABLE  V 


P 

P  c  -  c 

P  Psatd 

C 

P 

C 

P 

kg/cm^ 

atm  ^ 

cal/mole  “K 

Calc,  from 

HLJ  [35] 

Gutsche 
(smooth  curve 

7.7775  (V,P)  0 

11.74 

4o 

38.72  -5.63 

6.11 

7.05 

59 

57.1  -6.20 

5.54 

6.35 

79 

76.46  -6.51 

5.23 

5.70 

99 

95.8  -6.80 

4.94 

5.35 
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It  Is  interesting  to  note  that  Woolley,  et  al.,  (RP1932,  p.  U63)  con¬ 
cluded  that  Gutsche's  values  of  C  for  liquid  hydrogen  were  too  high 
to  he  consistent  with  the  best  thermal  and  P-V-T  data.  They  estimated 
that  Gutsche's  C  values  for  liquid  hydrogen  were  of  the  order  of  15 
per  cent  too  higS.  The  present  analysis  using  Johnston,  et  al.,  [35  1 
data  and  calcvilations  suggest  that  near  20“K  the  Gutsche  data  are  very 
neaurly  correct  but  do  Indeed  become  too  high  at  higher  teiq>erattire8. 

It  appears  that  a  plot  of  C  vs  T  over  the  range  20-50"K  and  P  = 
saturation  pressure  to  I50  Atmospheres  can  be  constructed  quite  readily. 
For  T  >  30“K  the  calculation  will  be  more  uncertain  largely  because 
C  for  the  liquid  varies  rapidly  with  temperature.  For  T  =  lU“  -  20" 
reliance  must  be  placed  on  the  limited  P-V-T  data  of  BarthdLonie  [lO  1  . 
Since  this  agrees  personably  well  with  the  data  of  Johnston,  et  ed.  [  9  ] 
at  20.38“K,  the  C  charts  could  probably  be  extended  to  -^lOO  atm  In  this 
region  as  well.  Values  of  C  could  also  be  calcvilated  from  the  NBS  T-S 
diagram  given  in  RP1932  by  the  relation 

0 

(If)  "  ^  isobar  (7) 

and  the  values  of  C  so  obtained  conqjared  with  those  calculated  from  the 
P-V-T  data  of  JohnsAon,  et  al. 

c.  Calculation  of  C  -  C  and  C 

p  V  V 


One  may  calcvilate  these  quantities  from  the  relations  given  In  RP1932 
as  follows: 


Cp  (p  ,T)  -  Cp"(T) 


=  -A-B  +  C-1 


(p  ,T)  -  C^»(T) 


=  -  A  -  B 


where  A  =  2  J  T  (||)^/p  3 

’■C 


Thus, 


C  -  C 

_E _ I 


=  C  -  1  + 


C  "  -  c  " 
P  V 

- R - 


=  C, 


C  =  C  -  RC 

V  p 


[Eq.  5.8  RPI932I 
[Eq.  5.7  RPI932I 


(8) 

(9) 


and 
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Chelton^  Maclnko  and  Dean  [^]  have  tabulated  the  quantities  necessary 
to  calc^ate  C  In  the  course  of  their  calculating  of  C  .  Their  values 
of  ^ 

’611  ,  been  used  to  obtain  by  equation  (9). 

The  final  results  are  shown  In  the  following  Table  VI. 
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TABLE  VI 

Ccnparleon  of  and  for  n-Hg 

P  =  4o  atm 


T 

"K 

P 

Amagat  density 

Cp 

Cal/mole*K 

C 

V 

70 

168.271 

6.545 

3.177 

60 

211.98 

7.426 

3.i4o 

50 

301.09 

9.751 

3.239 

k6 

368.81 

11.494 

3.223 

415.21 

12.106 

3.19B 

k2 

467.53 

11.754 

3.111 

30  atm 

60 

154.71 

6.832 

5.150 

295.96 

12.343 

3.282 

42 

346.96 

14.548 

5.196 

4o 

417.67 

15.366 

3.266 

38 

495.17 

14.672 

5.290 

20  atm 

44 

167.23 

8.935 

3.198 

4o 

219.29 

15.156 

3.322 

38 

278.49 

20.547 

5.475 

36 

418.30 

54.575 

3.58 

13  atm 

4o 

159.8 

9.013 

3.229 

36 

197.73 

16.852 

3.479 

3‘f 

410.66 

58.16 

3.18 
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These  calculations  shov  that  irtiereas  C  changes  rapidly  near  the  critical 
temperat\ire  for  pressures  less  than  UCratm.,  is  almost  constant  at 
3.2  cal /mole  *K,  very  nearly  the  same  value  as  C^". 

Direct  measiirements  of  C  for  gaseous  and  liquid  hydrogen  have  been  re¬ 
ported  by  Eucken  [26]  over  the  range  l8.6  -  47.5“K  for  a  few  cmstant 
values  of  the  Amagat  density  (60  -  86o).  !Rie  following  Table  VII  is  a 
summary  of  his  data  given  at  even  temperatures. 


TABLE  VII 

for  gaseous  cal/gfn  mole  "K 


T 

C  (moles /l.) 

2.67® 

4.89® 

16.95 

17.20 

22.30 

34.5 

36.2 

•K 

p  (Amagat 
density) 

59.9 

109.7 

380.1 

385.8 

500.2 

774 

612 

35 

3.20 

3.4o  ca 

4.5 

ca  4.0 

3.32 

3.39 

3.36 

37.5 

3.19 

3.35 

3.65 

5.50 

5.32 

3.42 

3.43 

iK) 

3.18 

3.28 

3.43 

3.39 

3.29 

^^5 

3.1^^ 

3.30 

3.30 

3.30 

3.68 

Earlier  work  by  Eucken,  Berl.  Sltzungsber.  (1912)  p.  l44. 

Eucken 's  results  are  for  the  most  part  In  very  good  agreement  with  the 
results  obtained  from  RP1932  and  given  in  the  preceding  Table  VI.  Eucken 
states  that  the  C  appears  to  pass  throu^  a  maximum  In  the  neighborhood 
of  the  critical  pSlnt  (”  in  agreement  with  researches  of  Dieter  id  on 
other  substances'). 


for  Liquid  Qydrogen 


C  has  been  measiired  for  liquid  hydrogen  by  Eucken  [26]  and  Bartholcne 
and  Eucken  [2^]  (according  to  RP19^)  for  a  series  of  densities  ranging 
from  380*^60  Ama^ts. 
expression 


The  relation  between  C  and  C  is  given  by  the 

p  V 


(10) 


Johnston,  et  al.,  [35]  have  tabulated  V,  ^  (^)  f 


-Mil) 


as  a  function  of  T  at  1°  intervals  from  20-33**^  and  for  pressures  from 
2>120  atm.  These  quantities,  together  with  the  C  results  obtained  earlier 
for  liquid  Hg  could  be  used  to  calculate  C^.  The^congnited  values  could 
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be  conpared  with  the  results  mentioned  above  (25#  26  1  . 

e.  Behavior  of  C  and  C  Near  the  Critical  Point 

_ £ _ 1 - 

[Dodge,  Chem.  Engineering  Iliermodynamlcs  p.  22Jff  gives  many  of  the  equations 
used  In  this  section]  . 


Uie  critical  point  is  ustially  characterized  by  the  two  relations 


Now  we  may  write 


(12) 


and  since  we  can  show  (see  below)  that  Is  finite  and^^j 

at  P  =  P  ,  T  =  T .  V  =  V  It  follows  that  (12)  that 
c  c  c 


(13) 


and  either  C  Is  finite  C  =  «o  or  C 
than  C^.  E^^rlmentally  Is  finite 


f. 


Uie  Nature  of 


is  an  Infinity  of  a  higher  order 
as  we  shall  see. 


P\ 

It  Is  an  experimental  fact  for  gases  that  Is  very  nearly  a  constant 

(l*e.,  the  Isometrics  are  linear)  for  all  V 


gases  examined.  The  usual  equations  of  state  reflect  this  fact.  For 
Instance,  for  the  van  der  Waal  equation  we  find 

R 

m 

■  V  -  b 

whence  f  d 

V 

=  0 

(15) 

For  the  Beattle-Brldgeman 

Equation  of  state  we  find  (Dodge,  p.  220,  eq.  ^) 

(S). 

V  ^ 

4  -  ^  (16) 

V 
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where  <>-  ,  P  and  X.  are  constants  at  fixed  T. 


o'  o 
We  see  that 


V  =  V 


We  may  also  note  that 


o 


=  finite  at  T  . 

c 


0  in  general. 


Frcn  the  relation 


gp' 


(17) 


we  conclude  that 


atT=T,P  =  P,V  =  V. 
c'  c'  c 


The  exact  behavior  of  C  for  z%al  gases  as  a  function  of  T  and  P  near 
the  critical  point  has  Seen  studied  experimentally  for  very  few  gases. 
The  classic  case  is  steam  [  See,  for  instance,  Dodge,  p.  226*7]  •  % 

has  been  examined  by  Chelton,  et  al,[  51  •  The  quantity  C  -  C  ° 
also  been  generalized  by  assuming  the  lav  of  correspondlng^stateS 
[See,  for  instance,  Hougen  and  Watson,  Chem.  Process  Principles,  Vol.  2, 
Dodge,  p.  2h3,  Lydersen  Greenkonaand  Hougen,  Eng.  Expt.  Station  Bull. 

No.  4,  Unlv.  of  Wisconsin  (1955) •] 


&  Behavior  of  in  the  Critical  Region 


We  may  write 


For  a  gas  obeying  van  der  Waal*s  equation  of  state  we  see  that  from 
Equation  0-5  ) 


=  •*,  T) 
T)  =  C 

for  all  T  and  V  Including  T  ,  V  . 

c  c 

Bence  =  f(T)  only. 


or 


(V,  T)  -  C°  (V  = 
Cy  (V, 


=  0 


o 


(V  =  ••  ,  T) 


(20) 
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For  the  Beattie-Brldgeman  equation  of  state  (Dodge,  p.  220)  c»ie  can  show 
that 


where  C,  B  ,  b  are  empirical  constants  in  the  Beattle*Brldgeman  equation 
of  state.  ^Consideration  of  the  magnitude  and  sign  of  C,  B^  and  b  for  some 
of  the  common  gases  [  see  Table  V-lt,  p.  18^  Dodge  ]  Indicates  that 

Is  small  and  negative 
Hence  C  -  C  °  >  0. 

V  V 

In  RPI952  one  finds  the  exact  expression,  expressed  In  terms  of  the  Amagat 
density 


C  (p,  T) 


^ - (  real  gas)  - 


C  °  (T) 

Y  ' 

R 


(ideal  gas)  =  -2 

(22) 


dp 

P 


Woolley,  et  al.  (  RP1932  )  have  tabulated  vs  p  and  T. 


We  note  that  J  >0  and  *•  ^  all  P  (0-500) . 

Hence  there  are  two  opposing  terms  In  Equation  (^2  ) . 


Writing  Eq.  (22  )  In  the  form 

*  lo  5  [-"  (P)p  -^0)^  1  P  ^ 

one  finds  that  A>0  for  0  ^  p  $  m  p  sr  55O)  A  =  0  «  p  s  350;  A <  0» 
for  350  ^  p  ^  500  for  T  =  36"K.  For  T  =  T0“,  200"K  .A>0  for  0  0  p  ^  500. 
[  These  sure  selected  T 's  for  calculation  1  . 


We  note  further  that  when  A  =  0 


T 


=  0,  In  this  case  a  maximum. 


(24) 
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The  consequence  of  this  is  that 


also. 


It  is  believed  that  a  detailed  analysis  of  C  vs  p  at  a  series  of  isotherms 
near  the  critical  (33*2U'*K)  will  shov  a  series  of  maxima  in  the  C  vs  p 
isotherms  for  Hg. 


Michels,  et  al.  [  32,  See  also  31 1  have  made  a  detailed  analysis  of  C  vs  p 
for  T  =  constant  for  C  of  COg  in  the  critical  region.  vas  calculated 
by  them  from  their  equation  or  state  data  for  COg.  T^se^cxirves  show  a 
maximum  In  C  for  COg  at  the  critical  point  (C  "  C  aril  at  T  and  p^). 
The  curious  ]Fact  pointed  out  by  Michels  is  that* the  maximum  occws  at  ^ 
very  nearly  the  critical  density  for  6lI1  Isothermt  frcn  t  -  31.0l»*  (=  t  )  to 
ten^eratures  as  high  as  130**C.  Michels  ccnnente  on  ihe  maxima  observed 
in  C  and  points  out  that  this  lollies 


for  p  ®  p  crlt. 


Michels  further  remarks  that  the  maximum  in  C  at  the  critical  density 
was  first  observed  by  Keesom  [32.1]  (this  azxlcle  is  not  available  to 
me)  and  remarks  that  Yound  observed 
2 

~  ^  Isopentane  of  the  critical  point .  [  32 .2  ] 


[  The  variation  of  C  with  p  at  constant  T  for  COg  is  also  discussed  in 
earlier  papers  by  Michels,  et  alJ  [  32.3^  32. U  ] 

It  is  interesting  also  to  note  that  Eucken  [26]  heis  observed  a  maximum  in 
C  near  the  critical  point  which  is  (as  he  remarks  without  reference) 

•Xn  agreement  with  the  researches  of  Dleterici  for  other  substances 
(Eucken' s  results  eu:«  given  in  Table  VII  of  this  paper.) 

h.  Generalization  of  C  -  C  ° 

V  V 


We  may  write  the  following  relations 

H  =  E+  PV  =  E  +  ZRT 

or  KPr 

c  c 

Since  the  terms  cm  the  rl^t  are  ejqjressible  in  terms  of  and  P_ 
(See,  for  Instance,  engineering  charts  such  as  given  by  Doiage, 
p.  241 ),  we  may  generalize 

E°  -  E 
T 

c 


(25) 

(26) 
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Finally  we  may  write 


'E°  -  E 


C  °  -  C 


'E°  -  El 


=  C  -  C 

V  V 


E°  -  E 


JS  *  £ 

Lydersen,  Greenkorn  and  Hougen  [  371  give  tables  and  one  plot  of - — 

However,  they  give  no  tables  or  plots  for  “  ^v* 

1 .  Determination  of  C  /C 

_  p  V 

Velocity  of  Soumd  Method.  Kie  velocity  of  sound  can  be  related  (for 
zero  frequency)  to  the  thermodynamic  properties  of  a  fluid  by  a  number 
of  relations,  such  as  (Zemansky,  p.  133) 


where  3  =*  - 


rnr 

V  IdP 


d  =  density  =  ^ 


We  may  also  write  for  C^/C  the  result 

p  V 


■  ^  ■  -kiU) 

s  -  1  /^V  x" 

V  Up  J 


These  equations,  together  with  other  thermodynamic  relations,  yield  the 
additional  relations 


=  -VIA  ( 

.2  -  .  jd  (n\  *  f.  JL. 

M  \9yJ^  M  c  VTt/, 
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The  relation  given  in  Equation  (50  )  is  Interest^g  in  that  it  relates 
the  equation  of  state  data  for  an  isotherm  to  a  ,  C  and  C  in  a  very 
clmple  way.  ^  ^ 


Equation  01  )  is  interesting  since  it  gives  a  simple  relation  between  § 

and  other  quantities  at  the  critical  point.  Thus,  when  T  =  T  ,  V  =  V 

c  c 


2 

^crlt 


are  finite. 


therefore 


(52) 


Resonance  Method  of  Measurement  of  C^/C  .  The  resonance  method  of 
_ p  V 

measuring  C  /C  is  described  briefly  in  Zemansky,  Heat  and  Thermodynamics, 
5rd  Ed.,  p.^l25ff.  It  apparently  does  not  require  any  P-V-T  data,  but 
is  a  direct  method.  Such  a  method,  if  feasible  for  low  temperature  use, 
woiild  be  very  useful  as  an  Independent  check  on  C  /C  or  as  a  means  of 
obtaining  C  from  measiired  C  data.  Several  references  to  this  method  are 
given  here  r27,  28,  29,  30  ]  T 


Thermodynamic  Properties  of  E^drogen 

Johnston,  et  al.  I  3^  have  calcxilated  the  quantities  C  -  C  H°  -  H  and 
S  -  S°  for  gaseous  hydrogen  from  their  Jovile-IHicmson  dSta  fBr  the  range 
65-300‘’K  and  for  pressures  from  2.5  to  200  atm.  Baehr  [21]  has  also  used 
the  Joule-Thomson  data  of  Johnston,  et  al.  [20]  to  calculate  C  -  C  ^  for 
the  temperature  range  60-85*’K  and  for  pressures  from  zero  to  ^  ^ 

l80  atm.  Baehr  has  also  calculated  the  adiabatic  and  Isothermal  J-T  co¬ 
efficients  over  the  same  ranges. 


The  heat  capacity  results  of  Johnston  and  Baehr  have  been  discussed  in 
the  section  of  this  paper  which  deals  with  heat  capacity. 

In  a  later  paper  Johnston,  et  al.  [3^]  have  computed  C  -  C  °,  H°  -  H  suid 
S  -  S°  from  their  very  extensive  P-V-T  which  extends  f^om  28-300®K  and  for 
pressures  of  one  to  200  atm.  The  results  are  presented  in  a  series  of 
charts  of  Z  vs  P,  T;  -  H  vs  P,  Tj  and  two  H  -  S  (MoUler)  diagrams. 
Unfortxmately  the  charts  are  difficult  to  read  accurately  because  of  the 
coarse  grid.  Tlie  C  -  C  results  have  been  discussed  in  the  section 
of  this  paper  whlch^deals  with  beat  capacity.  The  values  of  C  -  C 
given  by„ Johnston,  et  al.  are  certainly  Incorrect  in  the  range^*beloS 
about  42  K  for  pressures  above  40-60  atmospheres,  presumably  because 
the  virlal  equations  used  to  represent  the  P-V-T  data  in  this  region  do 
not  accurately  reproduce  this  data.  This  suggests  that  H  -  H°  and  S  -  S° 
may  also  be  in  error  in  this  same  region,  though  perhaps  to  a  lesser  extent 
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since  only  first  derivatives  are  involved  in  the  calculation  of  these 
quantities.  Since  the  P-V-T  data  of  Johnston,  et  al.  deviate  sigalficantly 
from  the  corresponding  values  of  Z  obtained  frcm  RP1952,  one  vlll  e}q>ect 
sane  disagreement  between  the thermodynamic  quantity  obtained  by  Woolley, 
et  al.  (RPI932)  and  Johnston,  et  al.  . 


Johnston,  et  al.  [351  have  also  used  their  P-V-T  data  for  liquid  n-H-  [9I 

J _ _ - _ J _ TTTT  ««  1  V  ^  ^ 


to  compute  V,  C  -  C 


H  -  H 


'  satd 


-  ^*2),  .  Igl) 


over  the  range  20-33"K  and  2-120  atm.  The  use  of  these  quantities  require 


that  C 


fsatei'  ^  'Batd 

is  discussed  in  the  section  on  heat  capacity.  A  detailed  comparison  of 


and  3 


satd 


be  established.  The  calculation  of  C 


these  quantities  with  RP1932  has  not  been  made. 
C  computed  with  the  help  of  Johnston's  C  -  C 


However  the  values  of 
values  agree  with 


*^satd 

experimental  error  with  the  measurements  of  Gutsche  l24l  at  20“K  for  the 
range  9-9^  atm.  When  RP1932  [l3  was  written  no  P-V-T  data  existed  for 
the  liquid  range  for  T>21‘’  except  cLLong  the  saturation  line.  It  appears 
that  the  data  and  computed  thermodynamic  functions  can  be  used  to  construct 
a  much  more  adequate  T-S  diagram  in  this  region. 
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